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ABSTRACT 

Context. Low-metallicity star-forming dwarf galaxies are prime targets to understand the chemical enrichment of the interstellar 
medium. The H i region contains the bulk of the mass in blue compact dwarfs, and it provides important constraints on the dispersal 
and mixing of heavy elements released by successive star-formation episodes. The metallicity of the H i region is also a critical pa- 
rameter to investigate the future star-formation history, as metals provide most of the gas cooling that will facilitate and sustain star 
formation. 

Aims. Our primary objective is to study the enrichment of the H i region and the interplay between star-formation history and metal- 
licity evolution. Our secondary objective is to constrain the spatial- and time-scales over which the H i and H n regions are enriched, 
and the mass range of stars responsible for the heavy element production. Finally, we aim to examine the gas heating and cooling 
mechanisms in the H i region. 

Methods. We observed the most metal-poor star- forming galaxy in the Local Universe, IZw 18, with the Cosmic Origin Spectrograph 
onboard Hubble. The abundances in the neutral gas are derived from far-ultraviolet absorption-lines (H i, C n, C n*, N t, O i, ...) and are 
compared to the abundances in the H it region. Models are constructed to calculate the ionization structure and the thermal processes. 
We investigate the gas cooling in the H i region through physical diagnostics drawn from the fine-structure level of C + . 

Results. We find that H i region abundances are lower by a factor of ~2 as compared to the H ii region. There is no differential deple- 
tion on dust between the H i and H ii region. Using sulfur as a metallicity tracer, we calculate a metallicity of 1 /46 Z G (vs. 1/31 Z 0 in 
the H ii region). From the study of the C/O, [O/Fe], and N/O abundance ratios, we propose that C, N, O, and Fe are mainly produced in 
massive stars. We argue that the H i envelope may contain pockets of pristine gas with a metallicity essentially null. Finally, we derive 
the physical conditions in the Hi region by investigating the Cn* absorption line. The cooling rate derived from Cn* is consistent 
with collisions with H° atoms in the diffuse neutral gas. We calculate the star-formation rate from the C n* cooling rate assuming that 
photoelectric effect on dust is the dominant gas heating mechanism. Our determination is in good agreement with the values in the 
literature if we assume a low dust-to-gas ratio (~2000 times lower than the Milky Way value). 

Key words, galaxies: abundances - HII regions - galaxies: individual: IZw 18 - galaxies: ISM - galaxies: star formation - 
galaxies: evolution 


1. Introduction 

Numerical simulations by Cen & Ostriker (1999), predicting 
the evolution of the metal content of the Universe, show that 
metallicity is a stronger function of density than age. It is there- 
fore plausible that gas-rich dwarf galaxies in the Local Universe 
may remain in an early stage of chemical enrichment, possi- 
bly preserving pockets of pristine gas (Kunth & Sargent 1986; 
Kunth et al. 1994). Blue compact dwarf (BCD) galaxies rep- 
resent a major testbed as they include the most-metal poor 
star-forming galaxies known (Kunth & Ostlin 2000). Among 
these, I Zw 1 8 (Mrk 1 1 6) continues to attract considerable in- 
terest and to feed intense debates in extragalactic research. 
For more than three decades, I Zw 1 8, discovered by Zwicky 
(1966) and first studied by Sargent & Searle (1970) held the 
record as the most metal-deficient galaxy known, with an oxygen 


* Appendix A is available in electronic form at 
http : //www . aanda . org 


abundance 12 + log (O/H) = 7.17 ± 0.01 1 in its NW component 
and 7.22±0.02 in its SE component (Izotov et al. 1999). Only re- 
cently has it been displaced by two more metal-deficient BCDs, 
SBS 0335-052W with 12+log (O/H) = 7.12 ± 0.03 (Izotov et al. 
2005) and DDO 68 with 12 + log (O/H) = 7. 14 ± 0.03 (Izotov & 
Thuan 2007). 

It was postulated that the difficulty of finding more metal- 
poor objects could be explained by thorough local mixing of 
the metals produced and released over a single starburst episode 
(Kunth & Sargent 1986). This hypothesis, implying the exis- 
tence of quasi-pristine gas beforehand, has been since chal- 
lenged by observations of uniform abundances in disconnected 
star-forming regions within a single dwarf galaxy (e.g., Skillman 
& Kennicutt 1993; Kobulnicky & Skillman 1997; Noeske et al. 
2000). Although evidence of local metal enrichment has already 
been observed (e.g., Walsh & Roy 1993), it is usually associated 


1 This translates into a metallicity of »1/31 Z 0 with the solar oxygen 
abundance taken from Asplund et al. (2009). 
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with the presence of Wolf-Rayet (WR) stars (Walsh & Roy 1989; 
Thuan et al. 1996; Lopez Sanchez et al. 2007). This hypothe- 
sis was later confirmed by Brinchmann et al. (2008), with the 
study of a large sample of WR galaxies from the Sloan Digital 
Sky Survey. 

Chemical abundances and metallicity are key parameters 
for the evolution of galaxies. Abundances in BCDs are usually 
probed in the ionized gas of the H n regions, as derived from op- 
tical emission-lines. While H n regions correspond to the present 
star-formation episode, they can also be related to previous 
episodes via the triggered star-formation mechanism, in which 
star formation take place following compression of molecular 
clouds by massive stellar winds and supernovae (e.g., Vanhala 
& Cameron 1998; Bhattal et al. 1998; Fukuda & Hanawa 2000; 
Pustilnik et al. 2001; Lee & Chen 2007; Elmegreen 2011). 
H it regions represent however a small fraction of the galaxy in 
mass as compared to the surrounding H i region (Thuan & Martin 
1981). 

The H i region provides another way of measuring the metal- 
licity, which accounts for the bulk of the mass of a galaxy. The 
H i region sets a reference abundance pattern, caused by metal 
enrichment over long time scales and large spatial scales, that 
needs to be understood and constrained. With this in mind, the 
comparison of abundances between the H i region and the H ii re- 
gions allows us to examine the dispersal and mixing of heavy 
elements from recent starburst episodes, and to constrain stel- 
lar evolution models. Furthermore, since future star formation 
thrives on the H i gas reservoir, the H t region makes it possible 
to study the interplay between star-formation history and metal- 
licity. While past and present star-formation episodes enrich the 
interstellar medium (ISM) in heavy elements, the latter are, in 
turn, able to facilitate and sustain star-formation through sev- 
eral gas cooling mechanisms (e.g., [Cn] 157 /mi, CO, and dust 
radiative cooling). The main questions that drive the studies of 
the Hi region in BCDs are the following. (1) What mechanisms 
control the Hi region metallicity? (2) What is the influence of 
metallicity on the star-formation history? 

BCDs are ideal targets because they display large amounts 
of Hi gas (Thuan & Martin 1981) and because the massive 
stars provide a strong far-ultraviolet (FUV) continuum on which 
absorption-lines from neutral species located along the line of 
sight are superimposed. Kunth et al. (1994) used the GHRS 
instrument on board the Hubble Space Telescope (HST) to 
observe interstellar absorption lines toward the many massive 
stars in IZw 18. The authors detected the Oi line at 1302.2 A 
and, together with the neutral hydrogen content measured from 
the 21 cm line, they estimated the metallicity of the neutral en- 
velope to be at least 10 times lower than that of the ionized 
gas in the Hii regions. Later, Thuan et al. (1997) obtained a 
GHRS spectrum for SBS 0335-052E and found the metallicity 
of the neutral gas to be similar to that of the ionized gas. These 
results remained however inconclusive because of possible satu- 
ration effects of the Oi absorption line (Pettini & Lipman 1995). 

A further step forward was then achieved with the launch 
of the Far Ultraviolet Spectroscopic Explorer (FUSE; Moos 
et al. 2000), which allowed the observation of Hi absorption- 
lines together with many metallic species such as N i, O i, Si n, 
P ii, Ar i, and Fe n. The spectral resolution of FUSE (A/AA ~ 
20 000) and its sensitivity made it possible to detect numerous 
absorption-lines, some of them apparently not saturated. The 
neutral gas chemical composition was derived in several BCDs, 
spanning a wide range in ionized gas metallicity, from 1/30 
to 1/3 Z Q , with IZw 18 (Aloisi et al. 2003, hereafter A03; 
Lecavelier des Etangs et al. 2004, hereafter L04), SBS 0335-052 


(Thuan et al. 2005), IZw 36 (Lebouteiller et al. 2004), Pox 36 
(Lebouteiller et al. 2009), Mark 59 (Thuan et al. 2002), NGC 625 
(Cannon et al. 2005), and NGC 1705 (Heckman et al. 2001; Lee 
& Skillman 2004). These investigations showed that the neu- 
tral gas of BCDs is not pristine; it has already been enriched 
with metals, with a minimal metallicity of — 1 / 50 Zq (see the 
FUSE sample analysis in Lebouteiller et al. 2009). The second 
most important result is that the metallicity of the neutral gas is 
generally lower than that of the ionized gas - except for the two 
lowest-metallicity BCDs, IZw 18 and SBS 0335-052E - imply- 
ing that the H i region has been probably less processed than the 
H n region. 

However the interpretation of these results has been limited 
by the lack of constraints on the properties of the gas detected 
in absorption in the FUV (physical location along the line of 
sight, density structure, temperature). FUV observations select 
out dust-free lines of sight toward the massive stars, and the 
intervening gas might have peculiar properties as compared to 
the global ISM. Uncertainties also remain in deriving column 
densities from absorption-lines in complex systems. Most FUSE 
studies of BCDs assumed a single absorption component aris- 
ing from a single homogeneous line of sight with the exception 
of Mark 59 (Thuan et al. 2002; see also Lebouteiller et al. 2006 
for the giant H ii region NGC 604 in M 33). A more realistic ap- 
proach would be to consider multiple lines of sight toward mas- 
sive stars, each line of sight intersecting clouds with possibly dif- 
ferent physical conditions (such as turbulent and radial velocity) 
and chemical properties. Whenever individual absorption com- 
ponents are unresolved, it is almost impossible to tell whether 
they are saturated, in which case the column density inferred 
from the global absorption line can be severely underestimated 
(see e.g., Lebouteiller et al. 2006). Following the suggestion of 
Kunth & Sargent (1986), Bowen et al. (2005) got around the 
problem by deriving the abundances in the neutral gas of the 
dwarf spiral galaxy SBS 1543+593 along a quasar line of sight. 
The authors found that the abundances agree with those in the 
ionized gas of the brightest H n region in the arms, contrasting 
with the FUSE results in BCDs. 

With the advent of the HST/Cosmic Origins Spectrograph 
(COS; Green et al. 201 1) we have the opportunity to look anew 
at the problem of the chemical discontinuity between the ion- 
ized and neutral gas in BCDs, with in particular the observation 
of the Sn multiplet, which Pettini & Lipman (1995) have pro- 
posed as a tracer of the neutral gas enrichment. The cornerstone 
galaxy I Zw 1 8 is therefore once again put to the test, not only be- 
cause of its extremely low metallicity content. As it turns out this 
object has led to discordant analysis of the same dataset obtained 
with FUSE by two independent groups of scientists. A03 used 
an individual line fitting method and found abundances globally 
lower in the neutral gas as compared to the ionized gas, while 
L04 find no significant difference between the two phases. L04 
find Oi/Hi = -4.7+°g, which is consistent with the O/H ratio 
observed in the Hu regions (-4.83 ± 0.03; Izotov et al. 1999), 
while A03 report a significantly different value with Oi/Hi = 
-5.4 ± 0.3, a discrepancy with the former result that awaits clar- 
ification. This is indeed worrisome since this object has also a 
considerable weight in the study of depletion and that of a pos- 
sible threshold for chemical enrichment. Most importantly, the 
discrepancy between the observational results complicates the 
interpretation of chemical evolution models. Recchi et al. (2004) 
found that abundances in the H l medium (defined by a tempera- 
ture lower than 7000 K) should be similar to those in the H ii gas. 
A definitive observational constraint is needed to challenge the 
models. 
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Fig. 1. HST/ACS image of IZw 18 (main body), with RBG colors de- 
fined as I, V , and B filters. Images were downloaded from the Hubble 
Legacy Archive (http : //hla . stsci . edu/) and later combined. The 
green contours show the H i column density distribution from Lelli et al. 
(2012) with 2" resolution (beam size in the bottom left). Contours are 
drawn for 3 (dashed), 6, 9, 12, and 15 x 10 21 cm -2 . The COS observa- 
tion, represented by the black circle (2.5" diameter), is centered on the 
NW star cluster. 

While providing accurate measurements of chemical abun- 
dances in the neutral gas, the COS observations also shed a new 
light on the physical conditions and thermal processes in the 
Hi region. The absorption line Cn* ,11335.7, arising from the 
fine-structure level of C + , is intimately related to the far-IR cool- 
ing line [Cn] 157 yum that was recently detected in IZw 18 with 
the Herschel Space Telescope (Lebouteiller et al., in prep.). We 
investigate in the present study the heating and cooling mech- 
anisms in the H i region in order to understand the present and 
future star-formation history of the galaxy. 

Our study is organized as follows. We present the observa- 
tions in Sect. 2. Column densities are determined in Sect. 3. After 
the description of the nebular models in Sect. 4, we examine the 
physical conditions in the Hi region in Sect. 5. We then derive 
and interpret elementary abundances in Sect. 6. Finally, we dis- 
cuss chemical evolution scenarios in Sect. 7. 

2. Observations and data processing 

2.1. Datasets 

To first approximation, IZw 18 is composed of two young, mas- 
sive star clusters, NW and SE, which ionize the surrounding 
medium. We used COS on the HST to observe the brightest 
region in the FUV, NW, in April 2010 (program 1 1 523) and 
later in December 2010 (program 12 028). Both observations 
were made with the circular aperture, with a diameter D ~ 2.5", 
corresponding to ~230pc at the distance of IZw 18 (~19Mpc; 
Aloisi et al. 2007). Observations were centered on NW (see 
Fig. 1 and the observation log in Table 1). Target acquisition con- 
sisted of a slew to the specified coordinates of NW followed by a 
peak-up in ultraviolet light. Although both observing programs 
made use of the FUV gratings G130M and G160M, the two 
sets of exposures cover somewhat different wavelength range. 
Together, they span ~ 1 135—1760 A. We also use the archival 
dataset from program 11 579 which observed with the G130M 


Table 1. COS datasets. 


# 

a 

(J2000) 

6 

(J2000) 

Grating//" 

11523 

09:34:02.00 

+55:14:28.0 

G130M/3.4, G160M/7.3, 
G185M/8.1 

12 028 

09:34:02.00 

+55:14:28.0 

G130M/8.2, G160M/1 1.4 

11579 

09:34:01.97 

+55:14:28.1 

G130M/15.5 


Notes. The exact coordinates for observations 11523 and 12 028 are 
unknown because a peak-up was performed beforehand to target the 
UV brightest region. Coordinates are J2000. (Q) Exposure time in units 
of xlO 3 s. 


grating (~1 132-1433 A). No peak-up was performed for the lat- 
ter observation and the telescope was slewed to the requested 
coordinates of NW. The G185M spectrum was not used because 
of the low data quality 2 . 

The pipeline reduction procedure is continually being im- 
proved, and we obtained our spectra well after the observations 
were made (July 2011) from the Hubble Data Archive, which 
makes an “on-the-fly” reduction using the best calibrations avail- 
able. We combined exposures from a given program to form a 
single, co-added spectrum. We did not merge the spectra from 
the three programs because the pointing was slightly different. 
This results in differences in the spectra, with for instance the 
Civ T1548, ,11550 doublet, which is formed in stellar winds and 
in the H u region. 

Among the 3 available datasets, the 1 1 579 observation has 
the longest exposure time below — 1433 A, while the 12 028 ob- 
servation has the longest exposure time above (Table 1 ). 
The 12 028 observations should be superior to the 11523 ob- 
servations by virtue of their longer exposure times. Indeed, the 
higher signal-to-noise ratio (S/N) in the 12 028 spectra enabled 
us to discern broad, shallow stellar absorption lines that we had 
not detected before in the other datasets. Also, the 12 028 ob- 
servations were accomplished in a single 5-orbit visit, so we 
are assured of a constant telescope pointing throughout the visit, 
whereas in program 1 1 523, the G130M exposures and G160M 
exposures were obtained in different visits. 

In some exposures, geocoronal emission (in particular from 
the Oi 1302, 1305, 1306 A multiplet) is contaminating the inter- 
stellar line profiles. We ignored the corresponding exposures for 
the wavelength range around the geocoronal emission. 

As a preparation for the fitting of the interstellar lines, we 
prepared each COS spectrum in several steps as illustrated in 
Fig. 2. In the four steps, we: 

- Corrected the observed spectrum for foreground Galactic ex- 
tinction by dust, using the color excess E(B - V) = 0.032 
(Schlegel et al. 1998). This correction is used for the nor- 
malization by the stellar model (see below); 

- Corrected for Lyo- absorption in the galaxy and for Ly a 
absorption by the high-velocity cloud in the line of sight 
to IZw 18, following Kunth et al. (1994); 

- Measured the column density of Ly a in I Zw 1 8 and removed 
Ly a (Sect. 3.4); 


2 The G185M spectrum contains in particular the Sin ,11808.0 line, 
with an oscillator strength about 50 times lower than the weakest 
Si ii line in our study, ,11546.7. However, the low S/N of the spectrum, 
due to the almost fully illuminated aperture (Sect. 3.2) and to cross- 
contamination between the specfi'al stripes of the grating on the detec- 
tor, does not allow setting a useful constraint on the Si n column density. 
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Fig. 2. Processing of the observed spectrum a) consists of correcting 
for reddening and for Lyu absorption in the Milky Way and in the high- 
velocity cloud b); fitting and removing Lya from IZw 18 c); and nor- 
malizing the spectrum d). The residual emission at Lyo- wavelength is 
due to geocoronal H i. 


- Normalized by the stellar spectrum (continuum and absorp- 
tion in photospheres; Sect. 2.2). Figure 3 shows the spectra 
along with the normalized stellar model. Final spectra are 
shown in the Appendix. 

2.2. Contamination by stellar absorption 

The FUV spectrum of I Zw 1 8 is clearly composite, with con- 
tributions from stars, H it region, and H i region. Absorption in 
stellar photospheres needs to be removed since they may con- 
taminate interstellar absorption-lines. We calculated model spec- 
tra of young stellar populations using grids of high-resolution 
(0.2 A) non-LTE photospheric spectra (900-2000 A) by Lanz & 
Hubeny (2003, 2007) and Geneva evolutionary tracks (Lejeune 
& Schaerer 2001). 

High-resolution imaging studies (Hunter & Thronson 1995; 
Heap et al. 2013) indicate that IZw 18 is undergoing continu- 
ous star formation rather than an instantaneous burst. Also, it 
is clear from deep HST/ACS observations that IZw 18 has a 
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widely extended, older stellar population (Aloisi et al. 2007; 
Contreras Ramos et al. 2011). However, it is equally clear 
that young, massive stars are the dominant contributors to the 
UV flux of the NW component. Observations of individual 
stars in IZw 18 with WFC2 (Hunter & Thronson 1995) and 
STIS (Brown et al. 2002; Heap et al. 2013) on HST indi- 
cate that the NW massive star cluster has a range of ages. 
The brightest stars in the UV are also the youngest, and their 
inferred masses range up to 150 M Q . We therefore assumed 
continuous star formation and a Salpeter initial mass function 
with lower and upper mass limits, 1-150 M Q , to calculate the 
spectrum of the stellar population for ages (i.e., the duration 
of star formation up until now) ranging from 0 to lOOMyr. 
We then convolved each calculated spectrum with a rotational 
broadening function (usin / = 150 kms -1 ) and with the instru- 
mental function (Sect. 3.2). From a comparison of these mod- 
els to photospheric lines in the COS spectrum, we derive a 
continuous-star-formation age of 10-15 Myr and a stellar metal- 
licity log(Z/Zo)~-1.7. 

Using the age and metallicity derived above, we calcu- 
lated the overall spectral energy distribution (SED) of stars 
in the NW component using the grids of low-resolution, 
LTE spectra by Castelli & Kurucz (2003). This grid and the 
TLUSTY/synspec grid only cover effective temperatures up 
to 50 000 K, somewhat lower than the highest temperatures of 
massive stars on the zero-age main sequence (65 000 K), and 
definitely lower than highly evolved stars such as WR stars. It 
is therefore likely that our model SED underestimates the ioniz- 
ing flux of the stellar population. 

Figure 3 shows the stellar models overplotted on each 
dataset. We find that the contamination from stellar lines is neg- 
ligible for most of the interstellar lines from I Zw 1 8 and also 
from the local absorption systems (Sect. 3.1). One exception is 
the contamination of S n 21250.6 in the Milky Way by the stellar 
C hi line. The stellar model was used to normalize the observed 
spectrum so as to isolate the purely interstellar absorption. 

3. Column density determination 

3.1. Method 

Once the spectrum has been corrected for stellar absorption 
(Sect. 2.2), several interstellar absorption line systems can be 
distinguished along the line of sight to the massive stars in 
IZw 18 (Table 2). The main interstellar lines are listed in Table 3. 
The atomic data, in particular the oscillator strength values, 
were taken from Morton (2003), except for P n and S n for 
which we used Federman et al. (2007) and Podobedova et al. 
(2009) respectively. The IZw 18 component is visible at a 
velocity of ~755kms -1 . Two local components are required 
at ~- 10 kms -1 and ~-90kms -1 . In addition, a component 
at »-60 km s -1 is detected that corresponds to ionized species in 
the Milky Way. Finally, the high-velocity cloud “A vi” in com- 
plex A is detected at ~ — 171 kms -1 (Hulsbosch & Wakker 1988; 
Wakker & van Woerden 1997). 

We noticed a relative shift in the velocity of the neutral 
component within IZw 18 between the datasets: 756 ± 3 kms -1 
for 12 028, 750 ±2 kms" 1 for 11 523, and 768 ± 3 kms” 1 
for 11579 (Table 4). The same relative shifts are observed for 
the Milky Way component, suggesting an instrumental origin. 
We attribute this difference to a different distribution of the UV 
flux in and outside the aperture. As a comparison to our velocity 
determinations, the Ha observation of the NW component, prob- 
ing the ionized gas in the H n region, is 742 ± 7 km s -1 (Petrosian 
et al. 1997) while the H i 21 cm line has a velocity of ~740kms -1 
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Fig. 3. COS spectra of IZw 18. The Ly a absorption (from IZw 18, the Milky Way, and the high-velocity cloud) was removed and the flux was 
normalized to the local continuum. The numerous absorption features often come in pairs: the weaker - redshifted - absorption feature is from 
IZw 18; the stronger feature arises in the Galaxy. A non-LTE, photospheric model spectrum (no winds) of the stellar population (lOMyr age with 
a continuous star formation hypothesis) and log (Z/Z 0 ) = -1.7 is shown in blue. While many features such as Chi ^11 175 and the Siiv 41400 
doublet appear to be well matched by the stellar model, other high-ionization features such as He n 41640 and the C iv 41549 doublet are not well 
matched. The origin of these lines (e.g., WR stars, Hu region...) is unclear. 


(Lelli et al. 2012). Table 4 also shows the radial velocities in- 
ferred from each species individually. 

The column densities were calculated with the line fitting 
routine Owens (Lemoine et al. 2002). Owens was originally de- 
veloped to solve the complexity of many interleaved absorption 
systems in the FUSE spectra. Although line blending is not a 
significant issue in the COS spectrum of IZw 18, we used Owens 


because of its ability to perform a simultaneous fit of the absorp- 
tion line profiles from any number of absorption systems. Each 
gaseous system is defined by (1) a radial velocity u rac i; (2) a tem- 
perature T; (3) a turbulent velocity b, and (4) a set of column 
densities from selected species. In practice, the species within 
a given absorption system are assumed to share the same set of 
parameters (u ral j, T, b). 


A16, page 5 of 32 


A&A 553, A16 (2013) 


Table 3. Main interstellar lines used in this study. 


Table 2. Absorption systems. 


System 

Species 

if 

(kms -1 ) 

b b 

(kms -1 ) 

IZw 18 

neutral' 

~755 

*17.8 

Milky Way 

neutral' 1 

*-10 

*22 


neutral' 

*-90 

*7 

Milky Way 

ionized (S m, Si m, P hi) 

*-60 

*85 

A vi cloud' 7 

neutral' 

*-171 

*21 


Notes. Ia> Radial velocity. Turbulent velocity. (c) The species in the 
neutral component are those listed in Table 5. (<7) High-velocity cloud. 


We assume in the following that the absorption arises from 
a single broad component as opposed to a collection of narrow 
components. Within this hypothesis, the column density deter- 
mined from saturated lines highly depends on the b parameter. 
Systematic uncertainties on b include errors on the instrumen- 
tal resolution (Sect. 3.2) or the fact that species in the neutral gas 
might not share exactly the same b value. If we choose to neglect 
these systematic uncertainties, the simultaneous fit of lines with 
a wide range of optical depths provides constraints on both b and 
the column density. We verified that the additional presence of 
saturated lines in the simultaneous fit does not bias the column 
density determination if the column density is otherwise already 
well constrained by weak lines. 

Other effects might also be at work that introduce system- 
atic errors on the column density. For example, the velocity 
distribution could include a number of individual narrow ab- 
sorption components that we cannot resolve. The integrated, 
unresolved, line profile may appear optically thin even though 
narrow components are saturated (the so-called "hidden satura- 
tion”). Although the absorption spectra of metal-poor objects, 
with their low metal column densities, is in principle less af- 
fected than more metal-rich sources, the hidden saturation effect 
is non-linear and can result in underestimating the total column 
density by several factors (e.g., Lebouteiller 2005; Lebouteiller 
et al. 2006). Observationally, this effect can be minimized by 
including lines with a low oscillator strength / (e.g., Cn*, Sii, 
P ii in our spectra). Another method is to verify that, for a given 
species, lines with a wide range of / values can be fitted by a 
single column density and that the apparent optical depth scales 
with / (e.g., N i, S n, Fe n in our spectra). Nevertheless, the col- 
umn density of some other species was determined only from 
strong lines. We note in particular Cn /I1334.5 (optical depth 
in the line center to ~ 76), Oi ,11302.2 (to ~ 34), and all the 
Sin lines, with the weakest being 41526.7 (to ~ 9). For these 
species, we rely on the hypothesis that the absorption does not 
arise from a collection of narrow components, and we keep in 
mind that the quoted errors could be underestimated. 

The absorption lines from the various systems are blended 
with each other throughout the COS range. The following 
lines in IZwl8 are contaminated by local absorption sys- 
tems, Fen 41142.4, Sin 41190.4, 41193.3, Mnn 41199.4, and 
Sii 41250.6. The Owens algorithm allows us to use these lines 
to determine column densities under certain conditions, either 
other lines from the same species are isolated and/or the column 
density of the local component is well constrained by other lines. 
In the FUV spectrum of I Zw 1 8, Fe n. Si u, Mn n, and S u have 
other isolated lines and the local absorption systems are well 
constrained by a variety of lines throughout the COS range. 

In order to derive column densities, Owens performs 
&X 2 minimization over a number of spectral windows containing 
one or several lines. While Owens enables the parametrization of 


Line 

4 

(A) 

/ 

Line 

4 

(A) 

/ 

Ci* 

1657.907 

4.71e-02 

Si 

1425.030 

1.25e-01 

Ci* 

1657.379 

3.56e-02 

Si* 

1323.516 

3.06e-02 

Ci 

1656.928 

1.49e-01 

Si 

1316.543 

3.45e-02 

Ci* 

1656.267 

5.88e-02 

Si* 

1310.194 

2.89e-02 

Ci* 

1560.682 

9.59e-03 

Si* 

1277.216 

4.93e-02 

Ci 

1560.309 

1.28e-02 

Si 

1270.780 

5.51e-02 

Ci* 

1329.123 

1.60e-02 

Si* 

1253.325 

2.90e-02 

Ci* 

1329.100 

2.60e-02 

Si 

1247.160 

3.24e-02 

Ci* 

1329.085 

2.13e-02 

S n 

1259.518 

1.66e-02 

Ci 

1328.833 

7.58e-02 

Sii 

1253.805 

1.21e-02 

Ci 

1280.135 

2.63e-02 

Sii 

1250.578 

5.43e-03 

Ci* 

1279.891 

1.43e-02 

Cl I 

1347.240 

1.53e-01 

Ci* 

1277.513 

2.10e-02 

Cii 

1335.726 

3.13e-02 

Ci* 

1277.283 

7.05e-02 

Cii 

1188.774 

7.01e-02 

Ci 

1277.245 

8.53e-02 

Mn ii 

1201.118 

1.21e-01 

Ci* 

1261.122 

1.47e-02 

Mnn 

1199.391 

1.69e-01 

Ci* 

1260.926 

1.35e-02 

Mnn 

1197.184 

2.17e-01 

Ci 

1260.735 

5.07e-02 

Mn ii 

1162.015 

1.02e-02 

Ci 

1193.030 

4.09e-02 

Fen 

1611.200 

1.38e-03 

Ci* 

1193.008 

2.91e-02 

Fen* 

1618.468 

2.20e-02 

Ci 

1157.910 

2.12e-02 

Fen 

1608.451 

5.77e-02 

Ci* 

1157.770 

1.83e-02 

Fen* 

1267.422 

1.69e-02 

Cii* 

1335.708 

1.15e-01 

Fe n* 

1266.677 

8.99e-03 

Cn* 

1335.663 

1.28e-02 

Fen 

1260.533 

2.40e-02 

Cii 

1334.532 

1.28e-01 

Fe ii’ 

1148.277 

8.28e-02 

Ni 

1200.710 

4.32e-02 

Fen 

1144.938 

8.30e-02 

Ni 

1200.223 

8.69e-02 

Fen 

1143.226 

1.92e-02 

Ni 

1199.550 

1.32e-01 

Fen 

1142.366 

4.01e-03 

Ni 

1134.980 

4.16e-02 

Fen 

1133.665 

4.72e-03 

Ni 

1134.415 

2.87e-02 

Niii 

1773.949 

6.22e-03 

Ni 

1134.165 

1.46e-02 

NiiC 

1754.813 

1.59e-02 

Oi* 

1304.858 

4.78e-02 

Niii 

1751.916 

2.77e-02 

Oi 

1355.598 

1.16e-06 

Nin 

1709.600 

6.88e-02 

Oi 

1302.168 

4.80e-02 

Nin* 

1788.491 

2.52e-02 

Ain 

1670.787 

1.74e+00 

Nin* 

1748.289 

1.40e-01 

Sin 

1526.707 

1.33e-01 

Nin 

1741.553 

4.27e-02 

Sin 

1260.422 

1.18e+00 

Nin 

1703.412 

6.00e-03 

Sin 

1193.290 

5.82e-01 

Ni ii 

1502.148 

1.33e-02 

Sin 

1190.416 

2.92e-01 

Nin 

1467.756 

9.90e-03 

Pi 

1787.648 

6.04e-02 

Nin* 

1500.434 

4.77e-02 

Pi 

1782.829 

1 . 13e-0 1 

Nin 

1467.259 

6.30e-03 

Pi 

1774.949 

1.69e-01 

Ni ii 

1454.842 

3.23e-02 

Pi 

1679.697 

5.36e-02 

Nin 

1415.720 

3.19e-03 

Pi 

1674.595 

3.25e-02 

Nin 

1412.866 

3.55e-03 

Pi 

1381.476 

3.16e-01 

Nin* 

1411.065 

6.28e-02 

Pi 

1379.428 

2.19e-01 

Ni ii 

1381.685 

5.08e-03 

Pi 

1377.073 

1.10e-01 

Nin* 

1381.286 

1.15e-01 

Pii 

1152.818 

2.72e-01 

Nin 

1370.132 

7.69e-02 

Si 

1473.994 

7.30e-02 

Nin 

1345.878 

7.69e-03 

Si* 

1433.309 

5.85e-02 

Nin* 

1335.201 

1.63e-01 

Si* 

1433.278 

1.69e-01 

Nin 

1317.217 

1.46e-01 

Si 

1425.188 

3.65e-02 

Nin 

1308.866 

5.41e-03 


Notes. 4 is the rest wavelength, / is the Uansition oscillator strength. 
The w denotes the transitions from the fine-structure level. 


the local continuum around the lines, the normalization we per- 
formed (Sect. 2.1) results in a flat continuum that does not need 
to be further parametrized. For each specie, we included every 
line in the COS wavelength domain. This includes the lines that 
are not detected, since the latter help providing a reliable upper 
limit to the column density. The fits were performed on the 3 ob- 
servation sets separately because of possible systematic differ- 
ences in the spectra. 
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Table 4. Radial velocities (2<x error bars). 


Species 

11523 

12028 

11579 

Global 

750 ±2 

756 ±3 

768 ±3 

Ci 




Ci* 




Cn 

755*5 

763!j 

769*? 

C n* 
Ni 

739 t 

735*3 

749*" 

737/5 

748!? 

169_ 5 

Oi 


753*5 


Al ii 

748*5 

765 !5 

Sin 

757 

763!f 

Sin* 




Pi 




Pii 




Si 




Sii 

747!" 

749+ 10 

763!f 

Cl i 




Mn ii 

749*5 

72i+ 10 

765 -6 

Fen 

759*? 

770** 

Ni ii 

733 " 2 

743!" 

769!" 


Notes. The first row gives the radial velocity of the neutral component, 
as derived from the simultaneous fit all the species in the table (as- 
suming that the species share the same radial and turbulent velocities; 
Sect. 3.1). For each species, we also provide the radial velocity calcu- 
lated by isolating the species in its own component, i.e., the radial ve- 
locity and the turbulent velocity are not constrained by other species. 

3.2. Line broadening 

The observed line width is a combination of the intrinsic broad- 
ening (thermal and turbulent) with the instrumental line spread 
function. The latter is a critical parameter upon which the b value 
determination relies. Systematic errors on the line spread func- 
tion could lead to wrong diagnostics on the line saturation and, 
therefore, to unreliable column density determinations. 

In the case of a point-source well-centered in the aper- 
ture, the instrumental full width at half maximum ( FWHM ) 
is ~0.08 A. If the source has extended emission or if it is a col- 
lection of point-sources, the spectral resolution is significantly 
degraded 3 . Depending on the source distribution in the aperture, 
the line profiles might be shifted and/or skewed, but we found 
the lines to be Gaussian and symmetric, suggesting a rather uni- 
formly distributed light in the aperture. 

We estimated the effective instrumental spectral resolution 
in two ways. First, we measured the one-dimensional flux pro- 
file of stars in the NW component as derived from a FUV im- 
age obtained by Brown et al. (2002). The profile looks roughly 
Gaussian with a FWHM = 1.25", i.e., half the size of the 
COS aperture. The flat-topped profile of geocoronal Ly a emis- 
sion, which uniformly fills the COS aperture, has a FWHM = 
0.84 A, so the expected instrumental profile is expected to have 
a FWHM = 0.4 A. The second estimate of the spectral reso- 
lution makes use of the data itself. For this, we considered the 
sharpest lines in the COS spectrum and also the line blends for 
which the 2 lines are barely separated (in particular at ~ 1134 A 
and ~ 1260 A). These profiles give a reliable constraint on the in- 
strumental FWHM , ~0.45-0.50 A. We use 0.50 A (~90kms -1 ) 
for the line fitting of the 11 523 and 12 028 datasets (Sect. 3.1). 
A slightly larger FWHM was required for the 1 1 579 dataset, 
0.55 A, which could be due to pointing jitter over the long 

3 COS instrument handbook, http://www.stsci.edu/hst/cos, 
Sect. 5.9. 


exposure time. For comparison, the point-like source FWHM 
with the FUSE telescope is 20 km s -1 (Hebrard et al. 2002; Wood 
et al. 2002). 

We used the instrumental resolution above to calculate the 
intrinsic line broadening. For the latter, we could not distin- 
guish between the thermal and turbulent components. We as- 
sumed that turbulence is the dominant source of broadening by 
setting the temperature to 1000 K, corresponding to a thermal 
broadening of <1 kms -1 . The turbulent velocity value we deter- 
mine for IZw 18 is b ~ 17.8 ± 3 kms -1 for the 12 028 dataset, 
and 18.2 ± 8 kms -1 for the 11 579 dataset, in good agreement 
with the value 17.3 ± 5 km s -1 found by L04. The 1 1 523 dataset 
did not converge on a well constrained b value, which is likely 
due to the relatively low S/N ratio of this observation and also 
to the presence of spurious features. We forced the b value in 
the 1 1 523 dataset to be ~ 17.8 ± 3 kms -1 . 

3.3. Heavy elements 

The fitting of the absorption line profiles is shown in the 
Appendix. Final column densities are given in Table 5. The error 
bars were calculated by using the Ay 1 method, which gives the 
variation of^ 2 -^A n , where ls l° west T 2 value corre- 
sponding to the global minimum over the entire parameter space. 
Since we use the information from all available lines in the wave- 
length range, the error bar can be significantly smaller than when 
considering a single line. We refer to L04 and Lebouteiller et al. 
(2006, 2009) for more details. The comparison between the col- 
umn densities we derive and the FUSE results asks for caution. 
The large FUSE aperture (30" x 30") includes the entire main 
body of I Zw 1 8 while COS only observed the NW component. 
However, most of the UV-bright stars are located toward NW 
(Brown et al. 2002), and the corresponding lines of sight pre- 
sumably dominated the FUSE spectral continuum. Therefore, 
we do not expect large differences between the column densities 
derived with COS and with FUSE. In the following we provide 
individual comments for each species. 

The C ii column density is not well constrained because the 
only Cn line, Al 334.4, is heavily saturated. Our determination 
relies on the confidence that the b value is well constrained 
by other lines from other species included in the fitting and on 
the fact that hidden saturation effects are negligible (Sect. 3.1). 
The column density of C n* is much better determined. C i and 
Ci* lines could be present in our spectra but they are either 
barely detected or not detected. The contribution from C i and 
C i* to the total C column density is negligible. 

The Ni ,11134 triplet, observed in the 1 1 579 dataset, is not 
saturated (see also L04), while the ,11200 triplet is close to satu- 
ration. The N i column density determination from the 3 datasets 
are in good agreement. They fall somewhat above the value from 
L04 and agree well with that from A03. 

Oi is observed in the ,11302.2 line (oscillator strength 
/ = 0.048) in the 1 1 523 and 1 1 579 datasets and in the A 1355. 6 
line (oscillator strength / = 10 -6 ) in all datasets. The former 
line is contaminated by geocoronal O i emission (Sect. 2.1), and, 
although the contaminated exposures were ignored, we cannot 
exclude that a weak contamination remains in the other expo- 
sures. The ,11355.6 line is never detected and provides a large 
upper limit to the Oi column density. L04 and A03 measured 
different O i column densities from the same spectrum, mostly 
because of the set of O i lines used for the fitting. L04 ignored 
the ,11039 line because of contamination by a terrestrial airglow 
and because of the line saturation. Our O i column density deter- 
mination is closer to that of A03, but we consider our value as 
highly uncertain because of saturation. 
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Table 5. Column densities. 




cos 


FUSE 

Model" 


11523 

12028 

11579 

Adopted 

L04 

A03 (±lcr) 


Hi 

21.34 

21.34 ±0.04 

21.34 

21.34 + 0.04 

21.30 ± 0.10 

21.35 ±0.10 

21.34 

Ci 


12-85!°$ 

<13.42 

<13.40 




Ci* 

<12.91 

<13.18 

<12.97 

<13.00 




Cn 

(15-46$$) 

Lf\ 
U> 
GC 
1 + 
© ° 
L> i© 

CTs OS 

(15.37$$) 

15.35 + 0.30 



15.49 (-0.05) 

Cn* 

13-84!°$ 

i n oq+ 0.14 
12).02)_oo7 

13-76$$ 

13.80 + 0.10 



13.53 (-0.37) 

Ni 

(14.55!°$) 

(14-45$$) 

(14.59$$) 

14.55 + 0.10 

14.22$$ 

14.44 ± 0.04 

14.54 (-0.01) 

Oi 

(16.05!°$) 

<18.08 

(16.00$$) 

-16.37* 

16.60$$ 

15.98 ±0.26 


Oi* 

(< 13.31) 


(<13.87) 

(< 1 3.90) 




Alii 

(13.32!°$) 

(13.25!°$) 

C 

13.32 ±0.30 



13.32 (-0.19) 

Sin 

(14.76!°$) 

(14-53$$) 

(14.51$$) 

14.52 ±0.30 

b b 
+ 

O 

OO 

14.81+0.07 

14.86 (-0.00) 

Sin* 

<12.31 

<11.69 

<12.06 

<11.90 




Pi 

<12.10 

<12.51 

<12.18 

<12.50 




Pii 

12 91 +0 - 21 
lz “ yi -0A9 

12.86$$ 

12 49 +0 - 31 

lz -^ y -l.03 

12.92 + 0.30 


<-13.60 


Si 

10 Q9+0.18 

1 ^’ yj -0.24 

<12.86 

<12.30 

<12.30 




S ii 

14 75+ 0 09 
1^+. / J _ 0 06 

4^ 

i + 
o © 
o o 

(14.56$$$ 

14.73 + 0.10 



14.74 (-0.00) 

Cl i 


<12.95 

<13.05 

<13.05 




Mnn 

13.10!°$ 

1 o 7T+O.I4 

<13.30 

13.10 + 0.50 



13.12 (-0.00) 

Fen 

(14.71!°$) 

(14-49$$) 

(14.76$$) 

14.75 ± 0.20 

14 60 +0 12 

i-+.uu _ 0 13 

15.09 ± 0.06 

14.66 (-0.00) 

Ni ii 

13.37!°$ 

1 O 4 72+0.08 

1 9 49+0-0 4 
-0.22 

13.42 ±0.10 



13.33 (-0.01) 


Notes. Error bars are 2<r error bars unless noted otherwise. Note that the error bars for the 1 1 523 dataset are likely underestimated because of 
the tight constraint on b for this dataset (Sect. 3.2). The values in parentheses could be affected by systematic uncertainties on the b value or due 
to spurious features (see text). The Oi column density in dataset 12028 is an upper limit since it is determined using the undetected 21355 line. 

Column density calculated by the Cloudspec model (Sect. 4). The value in parentheses gives the correction when the emission component is 
accounted for (Sect. 4.2). (i,) Estimated using S ii column density as a tracer (Sect. 6.2 ). (l?l No coverage . <d> Noisy profile. 


Aluminum has a single strong line in the COS range 
at 1670.8 A. The line is unfortunately saturated, and the 
Alii column density determination relies on the constraints on 
the b value. We adopt an average column density from the 
datasets 1 1 523 and 12 028. 

All the Sin lines are saturated. The Sin column determi- 
nation is thus relatively uncertain since it depends strongly on 
the b value determination and on the assumption that absorp- 
tion arises from a single broad component (Sect. 3.1). We use 
both 12 028 and 11 579 datasets for the Sin column density de- 
termination, and ignore the 1 1 523 dataset because of its rela- 
tively lower S/N. Our adopted column density is somewhat lower 
than the FUSE measurements by L04 and A03. We thus keep 
in mind that our Si n column density might be underestimated 
and adopt a large error bar that reconciles our value with the 
FUSE values. 

The strongest Pii line (21 152.8) is barely detected. Our col- 
umn density determination falls well below the upper limit from 
A03. The strongest Pi line at 1381.5 A is not detected, and the 
contribution from P i to the total P column density is negligible. 

Sii is observed in the 21250.8, 21253.8, and 21259.5 lines. 
The main constraint on the S n column density is given by the 
weakest line, 21253.8. It is well detected in the 12 028 dataset, 
while the profile is more noisy in the 2 other datasets, espe- 
cially in 1 1 579. We adopt an average column density between 
the 1 1 523 and 12 028 datasets. 

Chlorine is expected as Cli and Cl n in the neutral phase. 
Cl ii has no lines in the COS range and it was not detected in the 
FUSE spectra. Cli has numerous lines below ~ 1 400 A with the 
strongest at 1347.24 A. The latter is however not detected in any 
of the datasets. 


Manganese is observed via several Mn ii lines with strong os- 
cillator strengths. Unfortunately all the lines are heavily blended 
with saturated (or close to being saturated) N i lines from 
IZw 18 and from the Milky Way. Therefore, we consider that 
the Mn n column density is highly uncertain and may be affected 
by systematic uncertainties. 

Iron is observed via several Fe u lines spanning a wide range 
of oscillator strengths. The presence of weak lines should re- 
sult in a reliable Fen column density determination. Our col- 
umn density determination falls between the value from L04 
and A03. It agrees slightly better with L04 while it is not com- 
patible within errors with the value from A03. The difference 
between L04 and A03 is due to the lower b value (<8kms -1 ) 
assumed by A03 for the Fe n lines. 

Nickel is observed via many weak Ni ii lines, most of which 
are not detected. We stress that the column density determina- 
tion of Ni ii is better than 5cr because it uses the information 
from all the lines we included in the fitting. The column density 
determination from the 3 datasets are in remarkable agreement. 

3.4. Hydrogen 

The column density of H i is an important quantity since metal 
abundances are tied to it. However, it is a difficult parameter 
to measure correctly. The blue wing of the Ly a profile is con- 
taminated by numerous metal lines and by Lya absorption from 
the Galaxy and from a high-velocity cloud along the line of 
sight to IZw 18 (Sect. 3.1). Furthermore, the red wing of Lycr 
is contaminated by nebular continuum emission (2-photon and 
Balmer continuum) arising in the H n region of I Zw 1 8 and 
by the N v 21238, 21240 wind lines from the ionizing stars in 
IZw 18. Finally, the continuum level is uncertain, because Lya 
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Table 6. H i column density determinations. 


Value 

Method 

Instrument" 

Ref. 

Absorption 

21.34 ±0.04 

Lya (NW) 

HST/COS (=2.5") 

(1) 

=21.32 

Lyy3 

FUSE (30" x 30") 

(2) 

21.54 + 0.06 

Lya (NW) 

IUE & HST/GHRS (2") 

(3) 

21.34 + 0.07 

Lyman series 

FUSE (30" x 30") 

(4) 

21.32 + 0.09 

Lyman series 

FUSE (30" x 30") 

(5) 

22.3 

Lyor (peak) 

HST/STIS (52" x 0.5") 

(6) 

Emission 

=21.2 

2 1 cm (NW) 

VLA (8") 

(7) 

=21.20 

21 cm 

VLA (20") 

(8) 

21.53 + 0.15 

2 1 cm (NW) 

VLA (5") 

(8) 

=21.48 

2 1 cm (NW) 

VLA (2") 

(8) 


Notes. (a) The value in parentheses indicate either the aperture (absorp- 
tion measurements) or the beam size (emission measurements). 
References. (1) This study; (2) Vidal-Madjar et al. (2000); (3) Kunth 
et al. (1994); (4) A03; (5) L04; (6) Brown et al. (2002); (7) van Zee 
et al. (1998); (8) Lelli et al. (2012). 

is such a broad line. In principle, both red and blue wings can be 
contaminated by stellar Ly a absorption, but this contamination 
turns out to be insignificant (see also Lebouteiller et al. 2008). 

We worked with the COS spectrum processed through step 3 
(i.e., corrected for Lya in the Milky Way and in the high- 
velocity cloud; Sect. 2.1). We measured the Lya column density 
in IZw 18 by multiplying the observed spectrum by e T , with r = 
A(Hi) x 10 -19 - 37 x (A - A o)~ 2 (lenkins 1971). The column 
density that best canceled the observed absorption is N(H i) = 
2.2 x 10 21 cm -2 . We find a similar value for the 3 datasets. We 
list in Table 6 the various H i column density determinations in 
the literature, either derived in absorption from the Lyman series 
or in emission using the hyperfine-structure 21 cm line. Note that 
the comparison between observations in absorption and in emis- 
sion should account for the simple geometrical effect that the ab- 
sorption only probes the gas in the foreground of the UV-bright 
stars. 

The Hi spatial distribution observed by Lelli et al. (2012) is 
particularly informative. The high-spatial resolution (2" beam) 
reveals that the NW region lies in a Hi hole (see also Fig. 1). 
The NW region is located between the southeast H i peak, 
with N(Hi) ~ 1.5 x 10 22 cm -2 , and a region to the north 
with a rather flat column density distribution (~6 x 10 21 cm -2 ). 
Interpolating between these 2 values, we would expect an 
H i column density of ~10 22 cm -2 toward NW, which is signifi- 
cantly larger than our observed value. As discussed in Lelli et al. 
(2012), the young stellar cluster in NW likely created an H i hole, 
either due to ionization or to a blow out of the ISM. Assuming 
the H i hole is due to ionization, the column density of ionized 
hydrogen should be on first approximation 10 22 -2.2 x 10 21 = 
7.8 x 10 21 cm -2 . We discuss further the distribution of ionized 
gas in Sect. 4. 

4. Nebular models 

4.1. Geometry and physical conditions 

Models were constructed to compute the ionization structure of 
the nebula and the relative emission and absorption strength of 
the intervening medium along the line of sight to UV bright 
stars. We made the following simplifying assumptions: (1) the 
ionizing stars of NW can be treated as a single, central source; 
(2) nebular emission arises over the whole volume of H n region, 
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Fig. 4. Electron temperature (top) and density ( bottom ) profile of the 
simple model approximation of IZw 18. The abscissa is the radial dis- 
tance from the ionizing stars in cm. 

but nebular absorption arises only in material along the line of 
sight to the central source; and (3) there is a clear-cut boundary 
between the H n and H l regions. The latter assumption is prob- 
lematic. The photoionization model of IZw 18 by Pequignot 
(2008) indicates that the covering fraction of the Hu region is 
only 40%, so the bulk of the photoionizing radiation must es- 
cape into the H i region, and deep observations of I Zw 1 8 show 
Her emission in arcs outside the star-forming region and in a dif- 
fuse halo 3 kpc across (Dufour & Hester 1990). 

Spectral analysis of the Hi cloud of IZw 18 is much more 
complicated than for normal (i.e., star-less) damped Ly a sys- 
tems (DLAs). For example, the spectral profile of the C n 1334-5 
resonance doublet is formed by absorption by relatively lower- 
mass and/or cooler young stars, and by emission from the outer 
parts of the Hii region, as well as by absorption in the Hi cloud. 
Even Lycr, nominally formed in the H i cloud is not immune to 
contamination. Balmer continuum emission and 2-photon emis- 
sion from the H ii region affect the long-wavelength side of the 
line, and absorption by H° atoms in the atmospheres of cool stars 
can affect both the short- and long-wavelength wings of Ly a 
(Sect. 3.4). 

We used the software package, cloudspec (Hubeny et al. 
2000), which combines the photoionization code Cloudy 
(Ferland et al. 1998) with the spectral-synthesis code synspec 
(Hubeny & Lanz 1995) to compute the emission and absorp- 
tion spectrum of an object. The physical model is a medium 
photoionized by a central source and the cosmic background. 
Cloudy computes the ionization structure and emission spec- 
trum (both lines and continuum) of the medium. The synspec 
part of cloudspec then calculates the level populations and the 
attenuation spectrum produced by such a medium along the line 
of sight toward the central ionizing source. 

We ran cloudspec through the Hu region and into the 
H i region until the H i column density reached its observed value 
toward the NW component, 2.2 x 10 21 cm -2 (Sect. 3.4). Figure 4 
shows the radial profile of the electron temperature and density 
and the sharp change in the physical conditions between the H n 
and the H l region. Close in to the ionizing stars, the nebula is 
quite hot, due to the low abundance of oxygen, the main el- 
ement for collisional cooling in the ionized gas. The electron 
temperature decreases outward from the ionizing stars with an 
abrupt drop near the edge of the Hu region, and outside the 
Stromgren sphere, the temperature reaches down to 100 K. The 
electron density drops by 3 orders of magnitude at the edge of 
the Stromgren sphere, because most of the free electrons sup- 
plied by hydrogen are no longer available. 
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Fig. 5. Cloudspec model spectra. Top : stellar and nebular (H i and H u region) component normalized to the stellar continuum (blue) and emission 
component (red). Both component spectra have been convolved with a Gaussian line-spread function having a FWHM = 0.5 A to simulate the 
observations. Major absorption lines formed in the Hi region are labeled on top. Most other strong lines, e.g. the Siiv /11400 doublet, Civ .11549 
doublet, are higher ionization lines formed in the H ii region. The absorption component of Lyu completely erases the Lyu emission formed in the 
Hii region. Bottom : simple model spectrum for the case CF x AF = 0.1 (top) and for CF x AF = 0.4 (bottom), with CF the covering fraction 
and AF the aperture fraction (see text). The dotted line is the absorption spectrum, and the bold line is the net spectrum. The component and net 
spectra have been convolved with a Gaussian line-spread function having a FWHM = 0.5 A. 


The NW component is dominated by the H n region which 
sits an H i hole either due to hydrogen ionization or to a blow out 
of the ISM (Lelli et al. 2012). In our model, we find that the col- 
umn density of ionized hydrogen should be logATHi) = 21.79, 
i.e., almost a factor 3 larger than the observed H i column density. 
Following the hypothesis that the total (ionized+neutral) hydro- 
gen column density is flat across the main body of IZw 18, we 
would expect the ionized hydrogen column density toward NW 
to be around logiV(Hi) ~ 21.89 in order to compensate for the 
relatively lower Hi column density (Sect. 3.4). Our models are 
in good agreement with this value and support the idea that the 
gas was not blown out but mainly ionized. 

In assuming spherical symmetry, the models are highly sim- 
plistic. The ionizing stars in IZw 18 are offset from Ha emis- 
sion shell (Hunter & Thronson 1995; Brown et al. 2002), and the 


H n region has a knot of high density. Nevertheless, the models 
help us to understand important physical processes operating 
in IZw 18 (Sect. 5) and to guide the abundance analysis of the 
H l region (Sect. 6). 

4.2. Spectral line formation in the Hu region 

Singly ionized ions are formed by recombination of higher 
ionization stages inside the Hii region near the edge of the 
Stromgren sphere. Recombination and subsequent cascade down 
to the ground state in the H n region creates emission that works 
to reduce the apparent strength of absorption lines arising in the 
H i region. 

Figure 5a shows the emission and absorption components of 
the model line spectrum as derived by cloudspec. The observed 
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spectrum is the net. As can be seen from the figure, the emission 
is strong relative to the absorption for Cn* /11335.7. How im- 
portant is contamination by emission lines? The answer depends 
on two factors: the covering fraction of the Hu region, CF, and 
the fraction of the H n emission region subtended by the entrance 
aperture to the spectrograph, AF. Pequignot (2008) estimates the 
covering fraction of the H u region in I Zw 1 8 at CF = 0.4, based 
on the strength of H fi emission to Lyman continuum flux of 
the model best fitting the observations. The COS aperture frac- 
tion is not known, but is certainly much less than 1. Figure 5b 
compares a segment of the model UV spectrum for two cases, 
CF x AF = 0.1 (top) and CF x AF = 0.4 (bottom). It is clear 
that the latter case disagrees with observation by generally over- 
estimating the strength of emission, making the ground-state ab- 
sorption lines too weak, and the fine-structure lines appearing in 
emission (particularly C n* and Si ii*). The agreement is much 
better for the case where CF x AF = 0. 1 . In our abundance anal- 
ysis, we have therefore neglected the effects of emission except 
for the C ii doublet at 1334-5 A. 

Contamination by emission cannot always be neglected. In 
fact, the spectra of star-forming galaxies at high-redshift often 
show emission in the fine-structure lines (Pettini et al. 2000; 
Shapley et al. 2003). This emission is easily identifiable because 
these galaxies have outflows. The H i region is expanding with 
respect to the central, ionizing stars, so the absorption compo- 
nent, which arises along the line of sight to the stars, is blue- 
shifted. The emission component is formed over the outermost 
shell within the Hu region, so it is centered on the stellar ra- 
dial velocity. Expansion of the H u region simply broadens the 
line. Thus, the red half of the emission reaches the observer un- 
scathed, but the blue half of the emission is at least partially nul- 
lified by the blue-shifted absorption component. The net profile 
looks like a P Cygni profile. The emission spectrum of I Zw 1 8 is 
as real as in high-redshift, star-forming galaxies. The difference 
is that there is no detectable velocity gradient in IZw 18, so the 
emission is not separable from the absorption component. 

While the emission component in the H n region is negligi- 
ble for most lines, absorption can also arise from the Hn re- 
gion. A few resonant lines originating in the ionized gas are 
seen in absorption in the FUV spectra of I Zw 1 8 and provide 
an independent, direct, estimate of the ionization structure of the 
nebula. For instance, A03 measured the Fem T1122 line in the 
FUSE spectrum of IZw 18 and set an upper limit to the frac- 
tion of Fe ii in the ionized gas of 5%. A03 also cautioned that a 
large fraction of Sin absorption might arise in the ionized gas. 
We are now able to test this hypothesis, with the Sim A 1206.5 
line, which we detect at a velocity of ~755 ± 50 kins -1 in the 
COS spectra. This line partly originates from stellar atmospheres 
and partly from the Hu region. The observed profile of the 
Si hi line is complex and it depends on the stellar contribution 
that was removed (Sect. 2.2). The column density determina- 
tion also depends on the b value for which we have little in- 
formation in the ionized gas. Assuming that the line is not satu- 
rated, we find a consistent Si hi column density in the 3 datasets 
with ~2 x 1 0 1 3 cm -2 . What is the fraction of the Si n column den- 
sity we measured (3.3 x 10 14 cm -2 ) arising in the ionized gas? In 
the warm ionized medium in the Milky Way, about 90% of sil- 
icon is expected into Sin and 10% into Sim (Sembach et al. 
2000). Assuming this result is applicable to IZw 18, we find 
that only ~40% of the measured Si u column density might arise 
from the neutral gas. However, the situation is likely different in 
I Zw 1 8 where silicon in the ionized gas is expected mostly into 
higher ionization stages (Sect. 6.1; see also A03). We assume in 
the following that most of Si n should lie in the H i region gas. 


Unfortunately, no S m line could be used to perform a simi- 
lar analysis with sulfur. The brightest S hi line, A 1 190, is blended 
with Si n. Nevertheless, we consider that the ionization correc- 
tion smaller than for Si and that most of S ii lies in the H i region. 

5. Physical conditions and thermal equilibrium 
in the Hi region 

The C ii* column density represents only ~3% of the total car- 
bon column density (Table 5), yet Cn* is key to understand- 
ing the physical conditions, in particular the thermal processes 
toward the NW region. The ground state of C + is split into 
the fine-structure levels 2 Pi /2 and 2 P 3 / 2 - The absorption line 
Cn 41334.5 originates from the 2 Pi /2 level while Cn* T1335.7 
originates from -P3/2- The upper level is populated by colli- 
sions with H°, H 2 , and e~ . The spontaneous emission between 
the two fine-structure levels is observed in the far-IR emission 
line [C 11 ] 157 /.mi and dominates over collisional de-excitation 
as long as the density is lower than the transition critical density. 


5. 1. Electron fraction in the Hi region 


The population of the C + fine-structure level constrains the elec- 
tron density toward the NW region. Following lenkins et al. 
(2000a), the volume-density ratio n(C + *)/n(C + ) provides the 
electron density as long as collisions with electrons dominate 
the population of the fine-structure level. Collisions with H 2 
should be unimportant in IZw 18 given the lack of diffuse H 2 
(Vidal-Madjar et al. 2000). Furthermore, the collision strength 
with e - is larger than the collision strength with H atoms by 
a factor of ~150 for T = 6000 K. Collisions with e - will there- 
fore dominate as long as n e /n H remains larger than ~0.75% (e.g., 
Lehner et al. 2004). Potential sources of free electrons are soft 
X-rays, cosmic rays, and ionization of C°. We approximate the 
number density ratio n(C + *)/w(C + ) with the column density ra- 
tio Ren = n(C + *)/n(C + ) and, by assuming that electrons are the 
dominant collision partners, we derive an upper limit on the elec- 
tron density 


A 21 Ren 

721 (e) 2 e- r “c/r _ R cil 


-3 

cm , 


in 


where A 21 is the radiative decay probability from the fine- 
structure level from the 2 P 3/2 level to the ^Pl/ 2 > T eX c is 
the “P 3/2 level excitation temperature (~91 K), T is the electron 
temperature and 721 (e) is the collisional deexcitation rate coeffi- 
cient, which we take from Goldsmith et al. (2012). 

For a temperature range of [1000-8000] K, we find n e <, 
[10-20] Rciicm -3 - Using the observed column density ra- 
tio Ren = 2.8 x 10 -2 , we calculate n e < 0.4cm -3 . Since the 
C 11 absorption-line is saturated, we also checked that the result 
holds when using S 11 as a replacement for C n. The S 11 lines are 
not prone to saturation, and both S 11 and C n are the dominant 
ionization stages in the H 1 region, so we expect the ionization 
fractions of S into S + and C into C + to be well correlated over 
a wide range of conditions. Carbon is more readily depleted on 
dust grains than sulfur (Sofia et al. 1997, 1998), but rather than 
assuming a depletion factor, we assume that the abundance ratio 
in the neutral gas equates that in the ionized gas (C/S)hh ~ 0.98 
(Pequignot 2008). The latter assumption also implies that C and 
S disperse and mix in the ISM over similar timescales (Sect. 6.4). 
We can then approximate 7?cn with 


N( C + *) 1 

N( S + ) (C/S W 


( 2 ) 
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from which we derive a somewhat lower upper limit 
of « e ;£ 0.2 cm -3 . 

On the one hand, the upper limit on the electron density 
we determine is significantly lower than the value measured in 
the ionized gas of the star-forming regions derived in Dufour & 
Hester (1990) and in the Her arc derived by Izotov et al. (2001a), 
on the order of 10 cm -3 . This suggests that Cn* arises in the 
diffuse ionized gas or in the H l region. On the other hand, the 
upper limit on n e is significantly larger than the electron density 
in the Hi region inferred from our Cloudy models (~0.01 cm -3 ; 
Sect. 4). The models suggest that the Cn* emission component 
can be neglected and thus does not impact our estimate of 77cn 
(Sect. 4.2). However, it is plausible that collisions partners other 
than e - contribute to populating the C + fine-structure level, in 
particular H atoms. We can use the same method as Eq. (1) to 
derive an upper limit on the hydrogen density, i.e., assuming that 
hydrogen atoms are the dominant collision partners. For a tem- 
perature range of [100-1000] K, we find < [30-40] cm -3 . 

If we assume that the H i region has a density around or lower 
than 10 cm -3 , we must conclude that the electron fraction in 
the H i region is likely significant, possibly reaching a few per- 
cents. Such values would lie many orders of magnitude above 
the fraction of electrons provided by ionization of C° alone 
(equal to the carbon abundance relative to hydrogen), suggest- 
ing that the H i region is partially ionized, with electrons created 
by UV photons, soft X-rays, or cosmic rays originating from the 
star-forming regions. We investigate further the electron fraction 
in Sect. 5.2. 

5.2. Theoretical cooling rate from the Hi envelope 

In order to examine further the physical conditions in the H i re- 
gion, we now compare the observed cooling rate to the theoret- 
ical rate expected from the H i envelope, with H atoms and free 
electrons as collisions partners. The observed cooling rate per 
H atom is calculated following Pottasch et al. (1979) and Gry 
et al. (1992), with 

n( C + *) , n N( C + *) , , n 

L = A 2 i h V 2 i « 2.89 x 10 -20 — erg s -1 H -1 ) (3) 

wh (V(H u ) v ' 

where V 21 is the transition frequency and where we have again 
approximated the number density ratio with the column density 
ratio. We find Lqos ~ 8 x 10 -28 erg s -1 (H -1 ). We can rewrite the 
cooling rate per surface unit of the galaxy, 

L' ~ 2.89 x 10 -20 V( C + *) erg s -1 cm -2 , (4) 

which leads to L' cos ~ 1.8 x 10 -6 erg s -1 cm -2 . Assuming the 
source is extended, the intensity per surface unit of the galaxy 
is the same as the intensity per surface unit on the detector, so 
we can directly relate L' to the corresponding 157 pm emission 
intensity 7(157 pm) in erg s -1 cm -2 sr -1 via: 

7(157 pm) = H (5) 

4 n 

where <1> is the fraction of the beam filled with [Cn] emis- 
sion. Assuming <I> = 1, we find 7cos(157 pm) ~ 1.5 x 
10 -7 erg s -1 cm -2 sr -1 = 1.5 x 10 -10 Wm -2 sr -1 . This is the ob- 
served cooling rate toward the NW component. 

For the theoretical cooling rate in the H 1 envelope, we con- 
sider a medium with most free electrons produced by cos- 
mic rays, soft X-rays, or escaping UV photons (see Sect. 5.1). 


x= -4 -3 -2 -1 0 



n H [enr 3 ] 

Fig. 6. The 157 pm cooling line intensity is calculated from Eq. (6) and 
is plotted as a function of the hydrogen number density and of the elec- 
tron fraction ( x = n e /n H ). Note that the result should be scaled by the 
fraction of carbon in C + . The result should also be scaled by the carbon 
abundance with respect to the solar value and by the column density of 
hydrogen nuclei in units of 10 21 cm -2 . For temperatures 50 K and 100 K, 
only the curve for x = 0 is shown. 


Following Crawford et al. (1985) and Madden et al. (1993), the 
cooling rate expressed in W m -2 sr -1 can be written as 

7 th (157 pm) = 5.8 x 10 -7 f H F c+ ® y, (6) 

(Z Q ) (10- 1 cm - *-) 

where Z is the abundance of carbon relative to the solar value 
(C/H = 2.69 x 10 -4 ; Asplund et al. 2009), Nn is the hydrogen 
nuclei column density, Fq* is the fraction of carbon into C + ions 
along the line of sight, <f> is the fraction of the beam filled by 
C + gas, and y is the dimensionless factor defined as 

_ 2 e~ 9l ' T 

y 1+2 e -91/r + («e/«cr,e + « H /«cr,H) -1 ’ 

where n cr e and n cr ii are the critical densities for collisions with 
electrons and hydrogen (which we take from Goldsmith et al. 
2012). We assume in the following Fq* = 1 and we define x as 
the electron fraction x = n e /«H- Figure 6 provides 7 t h(157 yum) 
as a function of density for various temperatures and electron 
fractions. When the ionization fraction is large, proton impact 
could contribute to the 2 P 3/2 level excitation, but here we neglect 
this contribution since we focus on the H 1 region. It can be seen 
that 7* ( 1 57 pm) flattens for densities > 1 0 cm -3 for large x values 
(medium significantly ionized), and >10 3 cm -3 for low x values 
(purely neutral medium), which is a consequence of the criti- 
cal densities for collisions with electrons and hydrogen atoms. 
7 t h(157 pm) depends little on temperature for T greater than a 
few hundreds K. 

In IZw 18, the abundance of carbon relative to solar is Z = 
10 -2 , where we have used sulfur abundance as a substitute and 
assumed the C/S ratio in the ionized gas (see Sect. 5.1). We first 
consider a mostly neutral medium, with T ~ 1000K and use 
the H 1 column density derived in Sect. 3.4 (10 21,34 cm -2 ). From 
Fig. 6, we calculate that the observed cooling rate 7cos(157 pm) 
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can be reproduced by a wide range of conditions, from purely 
neutral with ~ 15 cm" 3 , to a medium with a significant elec- 
tron fraction with ~ 5 cm -3 and x ~ 2%. Assuming a ho- 
mogenous medium with a low clumping factor, such number 
densities imply small Hi scales on the order of >100 pc, which 
is compatible with an H t shell, left over after the NW star cluster 
blew out and ionized the pre-existing cloud (Sect. 4.2). 

For larger electron fractions, the medium is diffuse and ion- 
ized and we need to account for the column density of H + . 
Following the hypothesis that the total hydrogen column density 
toward NW is ~10 22 cm -2 (Sect. 3.4), we find that the cooling 
rate can be reproduced by an electron density of n t ~ 0.04 cm" 3 . 
However, such a low number density would result in an unreal- 
istic nebula size (>50kpc). 

In summary, although the lack of knowledge on the hydro- 
gen number density in the H i region prevents us from deriving 
a precise electron fraction, our results are compatible with most 
of the C ii* absorption arising in the H i region with likely a sig- 
nificant electron fraction. 

How does the cooling rate derived from Cn* compare 
with the predicted value from the Cloudy models? The mod- 
els (Sect. 4) predict that radiative de-excitation of the C + 
fine-structure level dominates the cooling in the Hi region, 
with -L mo dei ~ 3.5 x 10 -28 ergs -1 (H -1 ), which is a fac- 
tor of ~2 lower than the observed value Laos ~ 8 x 

10 -28 erg s -1 (H -1 ). While we assumed a constant hydrogen 
density of 10 cm -3 in the models, a lower density would re- 
sult in a higher temperature, with other coolants contributing 
significantly to the total cooling rate, such as [O i] 63 /rm or 
[Si ii] 35 yrm. This would suggest that the total cooling rate in 
the models is likely underestimated by several factors. Since our 
models assume a thermal equilibrium, the heating rate is also 
underestimated by the same amount. We discuss in Sect. 5.4 the 
heating mechanisms and how to reconcile the models with the 
observed C n* cooling rate. 

5.3. Cooling rate and SFR 

The observed cooling rate derived from Eq. (3) is Lqos ~ 8 x 
10 -28 erg s -1 (H -1 ). This value is significantly lower than in 
the diffuse medium of our galaxy (~10 -26 -10 -25 erg s -1 (H -1 ); 
Lehner et al. 2004) or in other nearby galaxies (e.g.. Madden 
et al. 1993), but it is close to the average cooling rate in DLAs 
(~10 -27 erg s -1 (H -1 ); Wolfe et al. 2003a). Assuming that the 
photoelectric effect on dust is the dominant gas heating mech- 
anism (see Sect. 5.4), the cooling rate scales with the dust-to- 
gas ratio (D/G), the photoelectric efficiency, and the FUV field. 
Since the latter is a good gauge of the star-formation rate (SFR) 
per unit physical area, the similar cooling rate between IZw 18 
and DLAs could in principle imply a similar surface SFR den- 
sity (Wolfe et al. 2003b). Following Bowen et al. (2005), we can 
estimate the SFR from C n* by using the Milky Way (MW) as a 
reference 

Ssfr(IZw 18) a [Z] -1 Xsfr(MW), (8) 

where Esfr is the SFR per unit area and [Z] represents the 
D/G normalized to the MW value. While [Z] should in prin- 
ciple scale with the gas metallicity, we use instead the D/G 
as measured with the Herschel Space Telescope (2000*}™ 
times lower than the MW value; Remy et al., in prep.). We 
assume Zsfr(MW) ~ 4 x 10 -3 M 0 yr -1 kpc -2 (Rana 1991) 
and L(MW) ~ 2 x 10 -26 erg s -1 (H) -1 (Lehner et al. 2004), and 


we finally obtain Zsfr(IZw 18) ~ 0.3 ± 0.2 M Q yr -1 kpc -2 . As 
an alternative measurement of the SFR, we use the empirical 
relation between the gas surface density S gas and Ssfr (e.g., 
Kennicutt 1998). The Hi column density is used as a tracer 
of Z gas , with a correction factor of ~4 due to the ionized gas 
along the line of sight (Sect. 4.1). We find X gas ~ 87.6 M 0 pc -2 
and Esfr(IZw 18)^ 0.13*}}'}}}} M Q yr -1 kpc -2 . Both these deter- 
minations are close to the values found by Aloisi et al. (1999), 
with 10 -2 - 10 -1 M 0 yr -1 kpc -2 , and Petrosian et al. (1997), 
with 0. 14 M 0 yr -1 kpc -2 (corrected to a distance of 19 Mpc). The 
slightly higher SFR found from the C n* cooling rate, if genuine, 
might be due to an underestimated D/G because of the presence 
of metal-free gas on the line of sight (Sect. 7.1) and/or to a larger 
photoelectric efficiency in I Zw 1 8 as compared to the MW (see 
also Sect. 5.4). 

5.4. Heating mechanisms 

In the warm ionized medium, photoionization of H should dom- 
inate the heating (see Pequignot 2008 for IZw 18), with likely 
extra heating from photoelectric effect on dust and from dissipa- 
tion of interstellar turbulence (Reynolds et al. 1999). The heat- 
ing mechanisms in the H i region are comparatively much less 
constrained. The gas heating rate in the ISM of our Galaxy is 
dominated by photoelectrons ejected from polycyclic aromatic 
hydrocarbons (PAHs) and small dust grains (Wolfire et al. 1995; 
Weingartner & Draine 2001), but low metallicity environments 
could show fundamental differences, with for instance the low 
PAH abundance of PAHs (e.g., Wu et al. 2006). 

Ultraviolet photons escaping the H ii region play an impor- 
tant role in the Hi gas heating through photoionization and 
through the photoelectric effect. Dufour & Hester (1990) discov- 
ered a faint halo of Her emission extending at least 2.7 kpc from 
the star-forming massive star clusters of IZw 18. Furthermore, 
Pequignot (2008) found that only 30-40% of ionizing photons 
from the NW and SE massive star clusters is absorbed by sur- 
rounding H ii regions. The leakage of ionizing photons into the 
neutral medium should result in a heating effect, affecting the 
relative population of the ground and fine-structure levels of C n, 
as well as the strength of the C n* T1335 fine-structure line. 

The total heating rate predicted by the Cloudy models 
in the Hi region (Sect. 4) is dominated by cosmic rays, and 
amounts to T ~ 3.5 x 10 -27 erg s -1 cm -3 . Photoelectric effect 
heating is comparatively negligible, with ~2% of the total heat- 
ing rate 4 . However, our simple model has important limitations. 
First, the heating rate is likely underestimated by several factors 
(Sect. 5.3), which could be due an underestimated photoelectric 
effect heating rate. The latter is highly dependent on the model 
geometry and parameters since it is proportional to the incident 
UV radiation field, the photoelectric efficiency, and the D/G ratio 
(e.g.. Bakes & Tielens 1994; Weingartner & Draine 2001). Apart 
from the remarkably low D/G in IZw 18 (Remy et al., in prep.), 
the low photoelectric heating rate predicted by our models is due 
to the relatively weak radiation field in the H i region (Go < 100, 
expressed in Habing units; Habing 1968), and to the low pho- 
toelectric efficiency, due to the large grain charging parameter 
(Go -\/7gas/we ~ 10 5 K 1/2 cm 3 ). The photoelectric effect might 
become the dominant heating mechanism in the models if a more 


4 If the D/G is simply scaled to the metallicity of IZw 18, the photo- 

electric effect heating rate becomes the dominant heating mechanism, 
contributing to half of the total heating rate, itself increased by a factor 
of 2. 
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realistic approach was taken. We list below some of the most im- 
portant parameters, 

- The radiation field in the Hi region depends on the 
ISM topology and in particular on the covering factor of the 
H ii region. The UV photons likely escape the star-forming 
region to reach and heat the gas on large spatial scales across 
the galaxy. 

- The electron density in the models might be underestimated 
(Sect. 5.1). A larger electron density will increase the photo- 
electric efficiency, resulting in a larger heating rate. 

- The grain size distribution could be biased toward the small- 
est sizes in low-metallicity dwarf galaxies (e.g., Galliano 
et al. 2005), with likely a significant impact in the gas heating 
rate and efficiency. 

- The D/G ratio is uncertain, between 700 and 3000 times 
lower than the MW value. 

The study of the far-IR [C n] line and its relation to the IR lumi- 
nosity will be important to confirm the importance of dust in the 
gas heating process (Lebouteiller et al., in prep.). 

Since our model appears to underestimate the heating rate 
(Sect. 5.3), we explored further the heating of the H i gas through 
the use of the “extra heating” parameter in Cloudy. The ex- 
tra heating represents all heating sources that are not well con- 
strained or not modeled (e.g., escaping ionizing photons, com- 
pressional heating, dissipation of turbulence, mechanical energy 
transfer from stellar winds and supernovae...). We assumed no 
radial dependence for this parameter. Figure 7 shows that the 
electron temperature, initially at 100K in absence of extra heat- 
ing, jumps from 300 K to 7000 K when the extra heating rate is 
increased from 10 -27 to 10 -25 erg s -1 cm -3 . The higher temper- 
ature brings increased collisional excitation of C + to the fine- 
structure level, and hence, stronger absorption of Cn* /I1335.7. 
As expected, the strengths of the emission components to the 
Cii doublet are hardly affected. The emission appears very 
strong - enough to overwhelm the absorption component of 
Cii /11334.5, but the COS aperture only subtended a small frac- 
tion of the Hii region, so the emission components are sig- 
nificantly weaker (see Sect. 4.2). Nevertheless, the unknown 
strength of the emission component adds to the uncertainty on 
the heating rate as derived from the observed C n* profile. 

Assuming an extra heating rate of ~10 -25 ' 8 erg s -1 cm -3 , we 
are able to reconcile our models with the cooling rate de- 
rived from Cii* (Sect. 5.3). The electron temperature T e in 
the Hi region with such extra heating increases from 100K 
to 1000 K. The electron temperature is unfortunately difficult to 
measure observationally, and, although the column density ra- 
tio N(Si + *)/N(C + *) provides in principle the gas kinetic temper- 
ature (Howk et al. 2005), our upper limit on the Si ii* column 
density is too large to provide any constraint. 

By including the extra heating, the main IR cooling lines pre- 
dicted by the models are [C it] (~40% of the total cooling rate), 
[O i] 63 jum (~30%), and [Si n] 35 pm (~26%). We wish to recall 
here that these values correspond to the diffuse medium, and that 
far-IR emission-line intensities might include other components 
(e.g., photodissociation regions). The Oi* 41304.9 absorption- 
line, arising from the 3 Pi fine-structure level of O 0 , is unfor- 
tunately undetected. Following Eq. (3), we calculate an upper 
limit on the [Oi] cooling of <10 _26 ergs _1 (H -1 ), which is sig- 
nificantly larger than the cooling by [C ii] inferred from C n* 
(Sect. 5.3). Therefore, we cannot examine in detail the effect of 
extra heating on the relative cooling line contributions. At best, 
the large upper limit on the cooling by [O i] is compatible with 
the extra heating. 
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Fig. 7. Effect of extra heating in the Hi region. The electron tempera- 
ture ( top panel) and the fraction of C ii and C n* to total carbon (middle 
panel) are shown as a function of radius to the ionizing stars. The pre- 
dicted line strength is shown in the bottom panel. The models assume 
unity for both the covering fraction of the H ii region and the fraction of 
the emission subtended by the COS aperture (see Sect. 4.2). 

6. Chemical abundances 

6. 1. Dominant ionization stages in the Hi region 

In order to calculate the abundances relative to atomic hydro- 
gen, all the ionization stages in the neutral phase should be ac- 
counted for. Fortunately, the COS range gives access to the lines 
of the dominant ionization stages in the neutral gas, namely C n, 
N i, O i, Mg it, Si ii, P n, S ii, Mn n, Fe n, and Ni n. Only chlo- 
rine is expected about equally between Cii and Cln. Most 
studies of DLAs assume that singly ionized metals such as C, 
Al, Si, S, Fe, and Ni reside exclusively in the Hi region so 
that elemental abundances can be obtained from the ratio of 
the ion column density to the Hi column density (e.g., Pettini 
et al. 2000). However, this assumption is not always correct for 
H i clouds with embedded star-forming regions. To illustrate, 
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we have constructed a simple Cloudy model of the gaseous 
medium surrounding the ionizing stars (Sect. 4). The model is 
simplistic in having a constant total hydrogen density, 10cm -3 , 
but in other respects - metallicity, ionizing radiation field, star- 
formation rate, and H i column density - the model is appropri- 
ate for IZw 18. 

Figure 8 compares the model spectral energy distribution of 
the ionizing stars incident on the surrounding gaseous medium 
and after passage through the Hu region. No H-, N-, or 
O-ionizing radiation gets through to the H i region, so oxygen 
and nitrogen are neutral in the H l region. However, stellar radi- 
ation longward of the H i Lyman limit, although diluted by dis- 
tance from ionizing stars, has basically the same relative flux 
distribution except for the addition of very strong Lyo- and other, 
much weaker nebular features. Hence, other elements such as C, 
Al, Ni, Mg, Fe, Si, S, and Cl will be ionized in the H i region. 

From Fig. 9, we can see that in the H i region, the dominant 
ionization state of metals with ionization edges longward of the 
H i Lyman limit is singly ionized. We use the column densities 
of the corresponding species to determine the total abundance. 
A fraction of these species also arises in the H n region where it 
can lead to emission (recombination of higher stages; Sect. 4.2) 
and absorption. 

A caveat should be that partially-ionized regions might ex- 
ist on the lines of sight that our model cannot account for. The 
existence of such regions might affect the ionization structure of 
species with a large photoionization cross-sections (PiCS). Even 
if the ionization potential is greater than 13.6 eV, the PiCS can 
be large enough for ionization to occur while hydrogen remains 
atomic. This effect is unimportant for thick clouds well shielded 
from the FUV radiation, but becomes significant for thin clouds 


or clouds embedded in a diffuse UV field. Among the species 
observed in absorption in galaxies, Ar i is the species the most 
affected because of its large PiCS (e.g., Sofia & Jenkins 1998; 
Jenkins 2013). More relevant to our study, N i and O i also have 
a PiCS larger than that of H i (Verner et al. 1996), but also lower 
than that of Ari by a factor of >3, so we expect a much smaller 
correction than for Ar. Most importantly, the ionization fraction 
of N and O are coupled to the ionization fraction of H via reso- 
nant charge exchange reactions (see also Jenkins et al. 2000b). 
For all the other species 5 , the first ionization potential is lower 
than 13.6 eV, and the second ionization potential is well above, 
so we always expect them as single-charged ions in the neutral 
phase even in partially-ionized regions. 

In the following, we compare our abundances to the solar 
values from Asplund et al. (2009) except for argon for which 
we choose Lodders (2003). For the Hu region abundances, we 
use the results of Pequignot (2008) and adopt error bars that en- 
compass the determinations by Izotov et al. (1999). Following 
Lebouteiller et al. (2004), we define the abundance difference 
between the H n and H i region as: 

S H ,(X) = log (X/H) h n - log (X/H) h „ (9) 

where (A/H) H n and (A/H) H i are the abundances of the element 
X in the H u region and H i region respectively. We also define 
the abundance relative to the solar value as: 

[A/T] = log(A/T)-log(A/y) 0 , (10) 

where ( X/Y) e is the solar abundance reference. 

5 Except Mn n, which is the only species for which we could not find 
the PiCS. 
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Figure 10 shows the abundances in the neutral gas (see 
Table 7 and Sect. 6) as compared those in the ionized gas. The 
line fits are shown in the appendix, including the line profiles as 
they would appear if the H t region abundances would be similar 
to those in the H it region. 

L04 and A03 examined the same FUSE spectrum and the dif- 
ferences are due to a different method of measuring the column 
densities and also to the choice of lines (Sect. 3.3). Our method 
is similar to L04’s method, i.e., a simultaneous fit of the lines. 
The differences in the abundances between our study and L04 
and A03 directly reflect the differences in the metal column den- 
sity determinations since our Hi column density is similar to 
their’s and since no ionization correction was used in our study 
nor in L04 and A03. 

6.2. O/H and the a-elements 

a-elements are produced in short-lived high-mass stars and are 
often used to trace the metallicity of galaxies. The majority of 
oxygen is produced by He and Ne burning during the hydrostatic 
burning phase with negligible contribution from the explosive 


stages. Silicon is produced by O burning (hydrostatic and ex- 
plosive) and Ne burning, while sulfur is produced by O burn- 
ing. Since the a-elements are made in the same massive stars 
that make oxygen, the abundance ratio between two a-elements 
should depend little on metallicity (see Izotov et al. 1999). 

The sulfur abundance is tightly correlated with that of oxy- 
gen in BCDs (Izotov et al. 1999). The S/O ratio in the ionized gas 
of IZw 18 is identical the solar ratio, with [S/O] = -0.07 ± 0.08. 
Within the error bars, we find that S is barely underabundant in 
the neutral phase with d Hr (S) ~ 0.18 ± 0.11. Figure A. 2 shows 
that the line profiles are not consistent with the H n region abun- 
dance. Sulfur is not expected to be depleted on dust (except 
in molecular clouds) and ionization corrections are negligible 
(Sect. 6.1), so we consider that this result implies a genuine un- 
derabundance of sulfur in the neutral gas. In other words, the 
metallicity of the neutral gas is between 51% and 85% that of 
the H ii region. 

The silicon abundance determination is fairly uncertain in 
the H ii region, with the values from Pequignot (2008) and Izotov 
et al. (1999) differing by almost a factor 2. Taking at face the 
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Table 7. Chemical abundances. 


X 

This study 

L04 (FUSE) 

A03 (FUSE) 

H ii region 

Adopted (X/H) Hn 

c 

6.02 + 0.30 



6.30 ± 0.10 (3), 6.55 (5), 6.54 + 0.10 ( 6 ), 6.47 + 0.15 (7) 

6.55 + 0.15 

N 

5.21+0.11 

4.92 + 0.20 

5.09 + 0.11 

5.58 ± 0.60 (1), 5.61 ± 0.03 (4), 5.60 (5) 

5.60 + 0.05 

0 

7.05 ± 0.15“ 

7 30+0.85 

6.63 + 0.28 

7.13 ± 0.05 (1), 7.17 + 0.03 (2), 7.22 + 0.01 (4), 7.21 (5) 

7.21 + 0.05 

Ne 




6.38 + 0.02 (1.2). 6.45 ± 0.02 (4), 6.39 (5) 

6.39 + 0.05 

Al 

3.98 + 0.30 





Si 

5.18 + 0.30 

5.50 + 0.30 

5.46 + 0.12 

5.61 + 0.23 (3), 5.89 (5), 5.65 ± 0.22 (7) 

5.89 + 0.30 

P 

3.58 + 0.30 


<4.25 



S 

5.39 + 0.11 



5.57 ± 0.40 (1), 5.75 + 0.06 (2), 5.59 + 0.05 (4), 5.57 (5) 

5.57 + 0.05 

Cl 

<3.71 





Ar 


4.15 + 0.30 

4.25 + 0.13 

4.87 ± 0.15 (1), 5.10 + 0.04 (2), 4.92 + 0.02 (4), 4.97 (5) 

4.97 + 0.10 

Mn 

3.76 + 0.50 





Fe 

5.41+0.20 

5.30 + 0.15 

5.74 + 0.12 

5.95 ± 0.90 (1), 5.41 ± 0.09 (4), 5.63 (5) 

5.63 + 0.20 

Ni 

4.08 + 0.11 






Notes. Abundances are given as 12 + log(X/H). w Calculated from Sn column density and assuming S/O = -1.64 (ionized gas value from 
Pequignot 2008). 

References. (1) (09 34 02.1 +55 14 25 J2000) Thuan & Izotov (2005); (2) (NW) Izotov et al. (1999); (3) (NW) Izotov & Thuan (1999); 
(4) (09 34 02.2 +55 14 25.06) Izotov et al. (1997), (5) Pequignot (2008), ( 6 ) Garnett et al. (1997), (7) Garnett et al. (1995). 



Ne Ar C N 0 S Cl P Mn Si Fe Ni Al 

Fig. 10. Abundance of each element in the H i region and in the H ii re- 
gion. Elements are ordered by increasing condensation temperature 
(temperature for which 50% of an element is removed from the gas 
phase). The horizontal lines show the metallicity for each phase based 
on fO/H], For the Hi region, S is used to estimate O (Sect. 6.2). 


value from Pequignot (2008), Si/O in the ionized gas is in good 
agreement with the solar ratio, with [Si/O] = -0.14 ± 0.33. 
Considering the similar nucleosynthesis process of these ele- 
ments, this indicates little or no depletion of Si on dust grains in 
the H it region. It must be noted that Peimbert & Peimbert (2010) 
used the relatively smaller silicon abundance found by Garnett 
et al. (1995) to conclude that ~39% of Si might be locked in 
dust in IZw 18. The silicon abundance in the Hi region is un- 
certain because of the Si n line saturation; our determination is 
somewhat lower than the column density derived by both L04 
and A03 (Table 5). Within errors, we find that silicon is under- 
abundant in the neutral gas with respect to the ionized gas, al- 
though with a large uncertainty, with <5 H i(Si) ~ 0.71 ± 0.42. 
Figures A. 2, A. 3 show that the line profiles are not consistent 
with the H n region abundance. Silicon can be affected by de- 
pletion on dust in the neutral phase, but this turns out to be a 
negligible effect (Sect. 6.7). Furthermore, although a fraction of 
Si ii could originate in the diffuse ionized gas (Sect. 4.2), it would 


Table 8 . Comparison between abundances in the H i and H n region. 


X 

(W 

[X/H] Hi 

[X/H] h „ 

6u,(X) 

c 

-3.57 ± 0.05 

-2.41 + 0.30 

- 1.88 ±0.16 

0.53 ±0.34 

N 

-4.17 + 0.05 

-2.62 + 0.12 

-2.23 ± 0.07 

0.39 ±0.14 

O 

-3.31+0.03 

-1.66 + 0.15" 

-1.48 ±0.06 

0.18 ± 0.16 

Ne 

-4.07 ±0.10 


-1.54 + 0.11 


Al 

-5.55 ± 0.03 

-2.47 ± 0.30 



Si 

-4.49 ± 0.03 

-2.33 + 0.30 

-1.62 + 0.30 

0.71 ±0.42 

P 

-6.59 ± 0.03 

-1.83 + 0.30 



S 

-4.88 ± 0.03 

-1.73 ±0.10 

-1.55 ±0.06 

0.18 + 0.11 

Cl 

-6.50 ±0.30 

<-1.79 



Ar 

-5.38 ±0.08 

-2.42 ± 0.17 c 

-1.65 ±0.13 

0.77 + 0.21 

Mn 

-6.57 ± 0.04 

-1.77 + 0.50 



Fe 

-4.50 ± 0.04 

-2.09 ± 0.20 

-1.85 ±0.20 

0.24 ± 0.28 

Ni 

-5.78 ± 0.04 

-2.14 + 0.11 




Notes. w Abundances in the present-day solar photosphere from 
Asplund et al. (2009) except for Ar for which we take Lodders (2003). 
® Calculated from S n column density and assuming S/O = -1.64 (ion- 
ized gas value from Pequignot 2008). (c) Argon abundance is taken from 
A03 and L04, with 12 + log(Ar/H) = 4.20 ±0.15 (Table 5). 

only increase the discrepancy between the silicon abundance in 
the two phases. The Si/S ratio in the ionized gas is solar ([Si/S] = 
-0.08 ± 0.30), while it is somewhat sub-solar in the neutral gas 
([Si/S] = -0.60 ± 0.32). Assuming that [Si/S] should also be »0 
in the neutral gas, we consider that our silicon abundance might 
be somewhat underestimated. This hypothesis seems to be con- 
firmed in Fig. 1 1, which shows that our observed [Si/O] lies be- 
low the predicted value from massive star yields (Woosley & 
Weaver 1995). 

The oxygen abundance is widely used as a metallicity tracer. 
However the O i lines in the COS spectra are either too weak to 
provide a useful upper limit or too strong to provide a reliable 
column density determination (Sect. 3.3). To circumvent this is- 
sue, we can use silicon or sulfur as a replacement. Sulfur is a 
better candidate because the S n lines are not prone to satura- 
tion. Bowen et al. (2005) find that Sii traces well Hi in a dif- 
fuse slab over a wide range of ionization parameter (see also 
Sect. 6.1). Sulfur is thus a reliable substitute for oxygen and 
a useful metallicity tracer. Assuming that the S/O ratio in the 
ionized phase of I Zw 1 8 (- 1 .64) also holds in the neutral phase, 
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Fig. 11. Abundance ratios calculated from the stellar yields of Woosley 
& Weaver (1995). Yields are plotted for two mass ranges, 11-25 M 0 
(open symbols) and 30-40 M 0 (filled symbols), and two metallicities, 
0.1 Z 0 (diamonds) and 0.01 Z 0 (squares). The grey stripe shows the 
observed values. The horizontal black line shows the average ratio ob- 
served in the H n region of BCDs (Izotov et al. 1999). 


we find (Os/H) ~ 7.03 ± 0.15. This value is barely in agree- 
ment with the O/H ratio measured in the ionized gas, leading to 
^Hi(Os) ~ 0.18 ± 0.16. Our indirect determination of the oxygen 
abundance lies between the FUSE value of L04, who found sim- 
ilar O/H between the two phases (although with large error bars), 
and the value of A03, who, despite using the same spectrum as 
L04, found b H i(0) ~ 0.54 (Sect. 3.3). 

6.3. The odd-Z elements 

Although Al and P are not a-elements, they are thought to be 
produced the same massive stars that make oxygen. No signif- 
icant amount should be produced in the explosive phases and 
the yields are expected to be metallicity dependent because of 
the odd-even effect (Woosley & Weaver 1995). Aluminum is 
produced by C and Ne burning while phosphorus is produced 
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during Ne burning. Lebouteiller et al. (2005) found that P n and 
O t column densities trace well each other over a wide range of 
environments as long as ionization corrections are not a major 
issue. Depletion on dust grains affects P more than O (Jenkins 
2009), so part of the dispersion observed by Lebouteiller et al. 
(2005) could also be due to a differential depletion factor. 

In this section, we discuss the abundances of P and Al rela- 
tive to a-elements since there is no measurement of P or Al abun- 
dance in the ionized gas. Assuming a solar P/O ratio we estimate 
(Op/H) ~ 6.86 ± 0.33. This value agrees within errors with the 
oxygen abundance (Os/H) ~ 7.05 ± 0.15 (Sect. 6.2). Figure 11 
shows that the expected yields from massive stars are also close 
to the solar P/O ratio. 

The Al/S ratio in the neutral gas is remarkably lower than the 
solar value, with [Al/S] = -0.74 ± 0.32. Aluminum abundance 
is determined using a single Al n line that is unfortunately sat- 
urated. However, the saturation of the Al n line is not as strong 
as C n or Si it lines, so we consider that saturation cannot explain 
the entire difference. It is also unlikely that there is a significant 
contribution from Al i given its extremely low ionization poten- 
tial of 5.99 eV. The low [Al/S] ratio could be due to a relatively 
stronger depletion of Al on dust grains as compared to S. The 
depletion of Al in the Local Interstellar Cloud is ~0.9 dex rela- 
tive to O in the WNM (Kimura et al. 2003). Such a factor would 
be more than enough to reconcile our Al/S ratio with the solar 
ratio. However, we investigate further the hypothesis of dust de- 
pletion in Sect. 6.7 and find that depletion alone cannot explain 
the low [Al/S] ratio. The Al/S ratio in the neutral gas, although 
sub-solar, is in fact consistent with expected yields from mas- 
sive stars (Fig. 11). We therefore conclude that aluminum in the 
H l region is not deficient with respect to the a-elements. 

6.4. Carbon 

Carbon is produced in intermediate- and high-mass stars by 
He burning. During the early stages of a galaxy evolution, car- 
bon is released mostly by the same young stars that make oxy- 
gen, so it is expected that C/O traces the yield ratio from massive 
stars. As galaxies evolve, older stellar population will contribute 
to the total carbon abundance, and a rise in C/O is expected (e.g., 
Izotov & Thuan 1999). 

In the ionized gas of IZw 18, [C/O] = -0.40 ± 0.17, which 
is somewhat lower than more metal-rich BCDs, suggesting that 
C mostly comes from massive stars, as proposed in Izotov & 
Thuan (1999). Note that the C/O measured by Izotov et al. 
(1999), -0.77 ± 0.10, is slightly lower than both the value 
in Garnett et al. (1997), -0.63 ± 0.10, and than the value in 
Pequignot (2008), ~-0.66. Garnett et al. (1997) found that C/O 
in IZw 18 is larger than expected from massive star nucleosyn- 
thesis and attribute the discrepancy to C enrichment by a lower 
mass stellar population. 

The carbon abundance in the H t region is somewhat lower 
than in the H n region. Using sulfur to trace oxygen in the neutral 
gas (Sect. 6.2), we find [C/Os] = -0.77 ±0.33, which falls some- 
what below, although compatible with, the ionized gas value. We 
also note that [C/Os] appears to be compatible with expected 
massive stars yields (Fig. 1 1). 

6.5. Nitrogen 

The nucleosynthesis of nitrogen is not as well understood as for 
carbon or oxygen. Nitrogen can be produced in the CNO cy- 
cle by any stars more massive than 1.5 M a . The observation of 
a N/O plateau at low-metallicity in BCDs (e.g., Izotov & Thuan 
1999; Nava et al. 2006) leads to the conclusion that nitrogen is 
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produced as a primary element, i.e., using the C and O produced 
in situ (since little C is available when the stars form). At higher 
metallicities (12 + log(0/H) > 8.3), secondary nitrogen pro- 
duction (from C and O from older generations) dominates. The 
origin of the N/O scatter at low-metallicity is still debated. If 
it is dominated by measurement errors, then nitrogen, like oxy- 
gen, is likely produced in short-lived massive stars (Nava et al. 
2006). The spread could also be intrinsic, due to the delayed 
injection of N from intermediate-mass stars. The observed scat- 
ter is theoretically explained either by effective yield variations 
and/or different star-formation history (Henry et al. 2000, 2006; 
Meynet & Maeder 2002b). The results of Henry et al. also show 
that nitrogen production in massive stars is possible. 

The N/O ratio in the ionized gas of IZw 18, -1.61 ± 0.06, 
is compatible with the low-metallicity BCD plateau located 
at ~-1.60 (Izotov & Thuan 1999). In the neutral gas, nitrogen 
is somewhat deficient, with <5 H i(N) ~ 0.39±0.1 1. This translates 
into a N/O ratio of N/Os = -1.84 ± 0.19 which is only some- 
what lower than the ionized gas ratio. Following the N/O plateau 
paradigm, we conclude that the ISM was enriched on large scales 
by the same stellar population. Assuming the lower N/O in the 
Hi region is genuine, it could be that intermediate-mass stars 
have contributed to the nitrogen enrichment in the H i region, or 
that nitrogen abundance was enhanced by massive stars in the 
H n region (see also Sect. 7.2). 

6.6. The iron group elements 

The iron group elements are particularly important because 
they could be released in the ISM on longer timescales than 
o'-elements. Iron is produced partly in massive stars and in 
SNela. A useful quantity to disentangle the two origins is 
[o'/Fe]. The [o-/Fe] ratio should increase with the progenitor 
mass (Woosley & Weaver 1995). If star formation is stopped, 
SNell contribute no more to the Fe enrichment and SNela 
could become important. Our Fe abundance is the same as the 
value determined by L04 using FUSE, and both values hint at a 
lower Fe abundance in the neutral gas as compared to the ion- 
ized gas. However, our value remains compatible with the ion- 
ized gas value within errors, with <5 H i(Fe) ~ 0.24 ± 0.28. We 
find that [Os/Fe] in the neutral gas is 0.45 ± 0.25, as compared 
to 0.37 ± 0.21 in the ionized gas. It is thus likely that the enrich- 
ment from SNe la, if any, is identical in the H i phase and in the 
H ii phase. This value of [Os/Fe] lies between 0 (solar O/Fe ra- 
tio) and the value in Galactic halo stars (Barbuy 1988) and in 
metal-poor stars in the solar neighborhood (lonsell et al. 2005), 
suggesting an early enrichment by high-mass stars (see Izotov & 
Thuan 1999). 

Nickel and manganese both belong to the iron group. The 
depletion strength of Ni is similar to that of Fe while Mn is less 
depleted than Fe by a factor <0.5 dex (Jenkins 2009). Manganese 
is mostly made in explosive Si burning while nickel is produced 
mostly by neutron capture on the iron-peak elements in SNe la. It 
should be noticed that Ni and Fe are among the most heavily de- 
pleted elements in the ISM (e.g., Savage & Sembach 1996). We 
discuss the influence of depletion on dust grains in Sect. 6.7. The 
Ni/Fe ratio in the neutral gas of IZw 18 is solar, with [Ni/Fe] = 
-0.05 ± 0.23. We note that [Ni/Fe] remains essentially null over 
a wide range of metallicities in nearby stars (also in DLAs; 
Izotov et al. 2001b) while at the same time [o'/Fe] increases at 
low-metallicity (Jonsell et al. 2005). This implies that Ni and Fe 
evolve in parallel, and it is therefore not surprising that we also 
find a solar ratio in IZw 18. Since we already suggested that 
iron might be produced mainly in massive stars, the same result 


should apply to nickel. Concerning manganese, its abundance is 
uncertain, and we measure a Mn/Fe ratio slightly above the solar 
ratio, with [Mn/Fe] = 0.32 ±0.51. 

6. 7. Depletion on dust grains 

We observe a small abundance discontinuity between the H i and 
the Hu region (Sect. 6). Could part of the discontinuity be ex- 
plained by a differential depletion factor on dust grains? Volatile 
elements such as sulfur and chlorine provide us with a way of 
estimating the fraction of the underabundance in the neutral gas 
that is not due to depletion on dust. Unfortunately, Cl is expected 
both as Cl i and Cl ii in the neutral phase so our upper limit on 
Cl abundance is not useful. On the other hand, sulfur is our best 
metallicity estimator since little or no S is expected to be locked 
in dust. We found that S/H is at least ~15% lower than in the ion- 
ized gas, implying that depletion alone cannot explain the lower 
metallicity of the neutral gas. 

We then examine whether depletion could explain the abun- 
dance discontinuity observed for the refractory elements. For 
this, we calculate the mass of each element in the Hi region 
locked in dust that would be required to explain the observed 
abundance discontinuity between the Hi and the Hii region. 
Given the small mass fraction of ionized gas in IZw 18 (e.g., 
Petrosian et al. 1997), we expect most of the dust mass to lie in 
the H i region and we therefore ignore the depletion in the H n re- 
gion. Dust extinction toward I Zw 1 8 is mostly Galactic, and the 
dust mass in IZw 18 is exceptionally low. Using data from the 
Spitzer Space Telescope, Galliano et al. (2008) found an upper 
limit of <2000 M 0 , while Herschel data constrains for the first 
time the dust mass, with 30CP™ M e (Remy et al., in prep.). We 
assume that the metals are well mixed in the H i region and ex- 
press the mass locked in dust m d (X) as: 

m d (X) = (/ - 1) m e (X) ~ (/ - 1) A(X) (X/H) obs M Hl , (1 1) 

where m d (X) is the mass locked into dust, m g (X) is the observed 
abundance in the neutral gas, / is the total mass (dust+gas) di- 
vided by the mass in the gas, A(X) is the atomic mass, (V/H) 0 bs 
is the observed abundance ratio, and Mm is the Hi mass. The 
total H i mass in IZw 18 is Mm = 6.9 x 10 7 M e , 38% of which 
(2.6 x 10 7 M 0 ) coming from the main Hi component associ- 
ated with the star-forming regions (Lequeux & Viallefond 1980; 
van Zee et al. 1998). 

The mass of each element in the gas phase is tabulated in 
Table 9. We also list the mass required in the dust phase to reach 
the corresponding abundance in the H ii region. For Al, P, and Cl, 
whose abundances are not constrained in the ionized gas, we es- 
timate the mass required to obtain a solar proportion with oxy- 
gen, i.e., [ X/0\ = 0. The mass of elements locked into the dust 
can be then compared to the dust mass measured by Herschel. 
Since we ignore the dust in the H ii region, the dust mass de- 
rived from Herschel is an upper limit to the available dust in the 
H l region. We find that depletion alone cannot explain the un- 
derabundance of C and Si in the neutral gas. Even though the 
required mass of N, Al, and S into dust is below the observed 
dust mass, altogether they represent at least half of the available 
dust, requiring unlikely dust compounds composition. We thus 
conclude that depletion on dust cannot explain the observed un- 
derabundances in the neutral gas. 

Since there is little dust available to remove elements from 
the gas phase, we also find that the abundance ratios are not 
affected by depletion on dust. As an illustration, there is not 
enough dust to reconcile O/Fe with the solar ratio (Sect. 6.6). 
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Table 9. Mass estimates of elements locked in dust. 


Element 

n h {xy 

Am* 


(M e ) 

(M 0 ) 

C 

325 +325 

-160 

357-776 

N 

59115 

59-86 

Os 

4670/}^ 

0-2100 

Al 

6 7 +67 

u - -3.3 

46-56 

Si 

110 /“° 

210-450 

P 

3+ 3 

-1.5 

0.4-3.4 

S 

2i4+ 61 

z,ah -48 

33-95 

Cl 

<4.7 

<4.4 

Fe 

074+220 
J 7 -138 

0-275 


Notes. (o) Mass of element in the gas phase of the Hi region (ob- 
served). (6) Mass required in dust to reach the Hii region abundance, 
i.e., <5 Hi (X) = 0. For Al, P, and Cl, for which there is no H ii region ref- 
erence, we provide the mass required to reach the metallicity calculated 
from O, [X/H] = -1.64. 

Similarly, we find that the low Al/S and Si/S ratio (Sects. 6.3, 
6.2) cannot be explained by aluminum and silicon locked in dust. 

7. Discussion 

The H t region of IZw 18 is not pristine; it has already been en- 
riched in metals. We also find that all the elemental abundances, 
in particular for Al, S, P, Fe, and Ni, are consistent with produc- 
tion in massive stars (Sect. 6). Our best tracer of the enrichment 
in the H i region by massive stars is sulfur (Sect. 6.2). By using 
the solar sulfur abundance from Asplund et al. (2009), we esti- 
mate the metallicity of the neutral phase to be « 1 /46 Z Q , as com- 
pared to ~1/31 Z;., for the H ii region. Depletion on dust grains 
should have a negligible impact in both phases considering the 
low dust mass in the galaxy (Sect. 6.7). We now discuss the en- 
richment of the H i and H ii regions. 

7. 1 . Abundance discontinuity H i-H ii 

The ionized gas in I Zw 1 8 is more metal-rich than the neu- 
tral gas, with abundances higher by ~0. 18-0.77 dex depend- 
ing on the element. Our results are a priori compatible with 
either (1) the Hu region being enriched by the present or re- 
cent star-formation episode(s) or (2) the presence of extremely 
low-metallicity clouds in the Hi region diluting the observed 
abundances. 

The elements synthesized and released by the current burst 
are expected to lie at first in a hot phase (10 6 K) and mixing 
timescales could be as long as 10 9 yr. Hot gas outflows have been 
observed through the O vi FUV doublet in starburst galaxies, in- 
cluding a few dwarf galaxies (e.g.. Grimes et al. 2009). However, 
there is no evidence that I Zw 1 8 developed such an outflow, even 
though the non-detection of the Ovi doublet could be due to 
the low oxygen abundance in this object (Heckman et al. 2002). 
Furthermore, Chandra X-ray observations of SBS 0335-052W, 
SBS 0335-052E, and I Zw 1 8 also failed to reveal any sign of hot 
gas breaking out from the stellar body implying that the galactic 
fountain scenario may not be valid in BCDs (Thuan et al. 2004; 
see also Sect. 7.3). 

As first hypothesized by Kunth & Sargent (1986), the Hu re- 
gion could be self-enriched, with a fraction of the newly pro- 
duced heavy elements dispersing on small spatial scales and 


mixing on short timescales (few Myrs). This would naturally ex- 
plain the H t-H ii region abundance discontinuity. Legrand et al. 
(2000) find that abundances are remarkably uniform in the ion- 
ized gas of IZw 18, with a dispersion of <0.1 dex on scales of 
several hundreds of parsecs, which seems to go against the fast 
and local enrichment hypothesis in the Hu regions. We recall 
that our H i region abundances are measured toward the NW re- 
gion, and compared to the Hu region abundances in that re- 
gion. The NW region is where the youngest stars are located, 
which in principle provides us with the best chance of observ- 
ing a local enrichment from the recent star-formation episode(s). 
The most significant abundance discontinuities H i-H ii we calcu- 
late are for N and Si, with <5 H i(N) = 0.39 ± 0.14 and h H i(Si) = 
0.71 ± 0.42 respectively (Sects. 6.5, 6.2), which is significantly 
larger than the abundance dispersion across the galaxy in its ion- 
ized gas. Although nitrogen abundance might be enhanced by 
very massive, rotating, stars in the Hii region (e.g., Meynet & 
Maeder 2002a; Hunter et al. 2008; Maeder & Meynet 2012), lo- 
cal enrichment by massive stars does not explain the abundance 
enhancement of the other elements in the H u region. 

Another possibility to explain a metal deficiency in the neu- 
tral gas is to invoke the presence of extremely metal-poor re- 
gions along the lines of sight, such as the low-metallicity high- 
velocity clouds infalling on the Milky Way. The H i envelope in 
IZw 18 extends as much as ~13kpc away from the star-forming 
regions (Lelli et al. 2012), and it is possible that the - presum- 
ably lower density - outskirts are not as enriched in heavy ele- 
ments as the immediate surroundings of the H ii regions. In order 
to estimate the mass of this relatively metal-poor H i gas, we ten- 
tatively assume that the H t medium can be separated into two 
components, one (hereafter Hl-sf) is spatially associated with 
the star-forming region and another component ( Hl-ex ) is more 
extended and possibly less enriched. In the following, we neglect 
the size and mass of the H it region. We further assume that the 
mixing of metals is uniform in the vicinity of the H n region, i.e., 
abundances are similar in the H ii region and in Hl-sf. The metal 
abundance in the H t region can then be written as 

1 + fx / H ex 

(X/H) Ul = / i ’ (X/H) sf . (12) 

1 + /H 

And the underabundance <5 hiTO in the neutral gas becomes: 

SufX) = - log [(1 - /h) + fx ./it, ex 1, ( 1 3) 

where / x = (X/H) ex /(X/H) s f is the relative abundance ratio of 
a given element between Hl-ex and Hl-sf and / H is the fraction 
of the total H i along the line of sight in the Hl-ex component. 
Figure 12 shows the variations of bmW with / H and f x . 

Taking sulfur as a metallicity tracer ( X = S), it can be seen 
that if Hl-ex is pristine (/s = 0), the observed sulfur deficiency 
is compatible with ~10M-5% of hydrogen into Hl-ex along the 
line of sight. Assuming a spherical distribution with constant hy- 
drogen volume density, this represents a factor [1 - ( 1 - /h) 3 ] in 
mass, translating to ~25-80% of the total H i mass in Hl-ex. In 
other words, about half the mass of the Hi envelope in IZw 18 
could be pristine. This scenario is therefore not compatible with 
a small number of low column density H i clouds such as the 
high-velocity clouds observed in the Milky Way. If the observed 
H i-H ii region abundance discontinuity is caused by the presence 
of pristine gas, the latter should represent a significant mass frac- 
tion of the galaxy (whose baryonic mass budget is dominated 
by H i gas). 

We can then assume that the Hl-ex component has also been 
somewhat enriched (/$ > 0). Considering that the H i mass in 
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f H 

Fig. 12. Metal deficiency in the H r region <5 Hi as a function of the frac- 
tion of hydrogen present in the Hl-ex component defined as the ex- 
tended H i component that may be less enriched in metals as compared 
to Hl-sf (the H i component associated to the star-forming regions). We 
assume that the Hl-sf component has the same metallicity as the H n re- 
gion and calculate <Shi for several values of the sulfur abundance ratio 
between Hl-ex and Hl-sf. The horizontal dashed line and the grey stripe 
indicates illustrate the sulfur deficiency, d'm(S) = 0.18 ± 0.1 1. 


Hl-ex is diluted over a larger volume, the aperture fraction is 
probably very small, and we can reasonably assume that the frac- 
tion of hydrogen in the Hl-ex component is less than 75%. For 
such a value of /fi, the observed sulfur deficiency is compati- 
ble with the sulfur abundance in Hl-ex being ~40-80% of the 
value in Hl-sf, and thus also of the value in the H n region. This 
would imply a metallicity as low as — 1/70 Zq in the outskirts 
of the galaxy. These results could put galaxy evolution scenar- 
ios to the test, with in particular the hypothesis of a minimal 
enrichment level at a given redshift (Sect. 7.3). Does a metallic- 
ity of — 1/70 Zq correspond to the minimal enrichment at red- 
shift z = 0? Deep Hi surveys such as ALFALFA (Giovanelli 
et al. 2005), and later SKA, should identify large numbers of 
optically-faint H i clouds in the Local Universe that will test the 
hypothesis of a minimal enrichment level at z = 0. We also wish 
to emphasize that the hypothetical presence of metal-free gas has 
important consequences on the interpretation of the integrated 
D/G ratio, in particular if the latter does not scale linearly with 
metallicity. 

7.2. Comparison with the sample of BCDs 

We now compare our results on IZw 18 to the sample of BCDs 
whose neutral gas has been investigated with the FUSE tele- 
scope. First, we examine the N/O abundance ratio as it poten- 
tially holds clues as to the stellar population responsible for en- 
riching the ISM. Nitrogen can be synthesized in a wider range 
of stellar masses than oxygen, and rapidly rotating massive stars 
can potentially enrich the ISM on short timescales (Sect. 6.5). In 
contrast, oxygen and the tr-elements are produced exclusively in 
massive stars and are released within a few Myrs through SNe 
(Sect. 6.2). Figure 13a compares the N/O ratio in the H i and in 
the H ii region of BCDs. The N/O ratio is similar in both phases 
for the highest metallicity objects. There might be a hint of a 
lower N/O in the Hi region for the lowest-metallicity objects 
(I Zw 1 8, SBS 0335-052, and I Zw 36), but the result is inconclu- 
sive within error bars. We conclude that the stars producing N 
have the same mass range in both phases, which argues for an 
enrichment of the ISM before the present starburst episode. 




- 2.5 - 2.0 - 1.5 - 1.0 - 0.5 0.0 

<[X/H]> hii 


Fig. 13. Top\ comparison of the N/O ratio between the H i and H ii re- 
gion for the sample of BCDs studied with FUSE, and for IZw 18 
with COS. The labels are defined as la) IZw 18 (L04), lb) IZw 18 
(A03), lc) IZw 18 (this study, with S as a tracer of O; Sect. 6.2), 
2) NGC 1705 (Heckman et al. 2001), 3) SBS 0335-052 (Thuan et al. 
2005), 4) Mark 59 (Thuan et al. 2002), 5) NGC 625 (Cannon et al. 
2005), 6) IZw 36 (Lebouteiller et al. 2004), 7) Pox 36 (Lebouteiller 
et al. 2009), 8) NGC 604 (Lebouteiller et al. 2006), 9) SBS 1543+593 
(Bowen et al. 2005; Schulte Ladbeck et al. 2004). Bottom', ratio of the 
average metal abundance [M/H] = ([X/H]) between the Hi and Hii re- 
gions. [M/H] is calculated two ways, averaging all elements observed in 
each phase (small diamonds), and averaging only the elements observed 
in both phases (large diamonds). 

The hypothesis of relatively metal-poor outskirts (Sect. 7.1) 
can also be put to the test in the BCD sample. Lebouteiller 
et al. (2009) summarized the FUSE results and showed that 
the metallicity in the Hi region in BCDs was equal (IZw 18, 
SBS 0335-052) or lower (IZw 36, Mark 59, Pox 36, NGC 625, 
NGC 1705) than that in the Hu region. A trend might exist, 
with the discontinuity being essentially null for the lowest metal- 
licity objects, and then constant around ~ 1 dex for metallici- 
ties >1/6 Zq (value in the ionized gas). However, if nitrogen 
is used as a metallicity tracer instead of oxygen, the results 
could also be interpreted as a constant discontinuity indepen- 
dent on metallicity. The lowest-metallicity object, IZw 18, is 
key to conclude on the ubiquity of the abundance discontinu- 
ity H i-H ii region. While the result for I Zw 1 8 was at the time 
inconclusive because of the discrepancy between A03 and L04, 
our study suggests that the abundance discontinuity in I Zw 1 8 
is rather small, smaller than what is observed in other BCDs. 
We now revisit the findings of Lebouteiller et al. (2009) by in- 
cluding this new measurement, and by calculating the average 
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metal abundance \M/ H] = ([X/H]) in the Hi and Hu regions 
(Fig. 13b). The choice of averaging the abundances is motivated 
by (1) possible systematic problems in the abundance determina- 
tion of some elements (hidden saturation, stellar contamination, 
atomic data...); and (2) the fact that there is so far no clear ev- 
idence of a different nucleosynthesis origin of elements in both 
phases. Figure 13b shows that the abundance discontinuity is 
around 0.7 dex for most objects, with no clear correlation as a 
function of metallicity. This implies that the chemical enrich- 
ment H i and H u region is disconnected in most objects, with the 
Hi region systematically more metal-poor than the Hu region. 
Applying the results of Sect. 7.1 to the sample of BCDs, Fig. 12 
shows that for hm ~ 0.7 dex, the Hl-ex component (extended H i 
envelope) would have a metallicity between 0% and 30% that of 
the Hl-sf component (associated to the H n region), and would 
represent more than 75% of the Hi mass. This result is rele- 
vant in particular for SBS 0335-052, which, despite a metallicity 
close to IZw 18, shows a discontinuity of ~0.7 dex (Fig. 13b). 
For SBS 0335-052, the Hi envelope could be as metal-poor 
as ~ 1/150 Z Q . New abundance determinations in this object by 
COS should bring a definitive answer (James et al., in prep.). 

7.3. Enrichment of the H i region 
7.3.1. Cosmic evolution 

Despite the weak discontinuity between the abundances of the 
Hi and Hu regions of IZw 18. the ISM is fairly well mixed 
on kiloparsec scales, with the metallicity of the neutral gas at 
least ~50% of that in the ionized gas (most likely value of 66%). 
How was the ISM, and in particular the Hi region, enriched? 
Considering the relatively uniform abundances, at least half of 
the metallicity in the H i region was presumably set before the 
current starburst episode. A major uncertainty concerns the gen- 
uine age of the galaxy as compared to its chemical age. The latter 
is probed by the metallicity parameter, which is determined by 
a complex interplay between successive star-formation episodes 
and external gas flow exchanges (outflow of enriched material 
and infall of metal-poor gas). 

Although the age of the oldest stars in IZw 18 is still sub- 
ject to debate (see Papaderos & Ostlin 2012), stars as old as 
the Hubble time (Aloisi et al. 2007; Contreras Ramos et al. 
2011) might be present, suggesting that IZw 18 is not genuinely 
young, and therefore that its present-day metallicity is the result 
of the mechanisms mentioned above over cosmic times. For this 
reason, the comparison between IZw 18 and high-;; objects is 
relevant and instructive. The intergalactic medium (IGM) was 
already enriched to metallicities similar to that of the Hi re- 
gion of IZw 18 at z ~ 1.6-2. 9 (Teller et al. 2002). Hence it is 
possible that IZw 18 was assembled from a collection of low 
column-density clouds in the IGM a few Gyrs after the Big 
Bang, and experienced little enrichment since then. The metal- 
licity of the Hi envelope of IZw 18 is also typical of DLAs 
at z > 2 (Rafelski et al. 2012), again compatible with a slow 
chemical enrichment. The lack of strong constraints on the mass 
of DLAs makes the comparison with I Zw 1 8 relatively uncer- 
tain within the chemical downsizing scenario in which galax- 
ies with different mass have different star-formation efficiencies 
(e.g.. Brooks et al. 2007). Nevertheless, the lack of nearby galax- 
ies as metal-poor as high-z DLAs (~ 1/100 Z Q ; e.g., Prochaska 
& Wolfe 1999; Rafelski et al. 2012) implies that nearby BCDs, 
including I Zw 1 8, have experienced a slow, minimal, enrich- 
ment over several Gyr. Legrand et al. (2000) suggested that 
I Zw 1 8 went through a mild continuous star-formation rate 


( SFR ), ~10“ 4 M Q yr -1 , over several Gyrs as opposed to a se- 
ries of violent short bursts separated by relatively longer quies- 
cent eras. Metal enrichment in BCDs through a series of violent 
star-formation episodes is an unlikely process as galaxies with 
different star-formation histories would result in different metal- 
licities, thus including a class of nearby galaxies with no star- 
burst episodes and an extremely low-metallicity, which is not 
confirmed by observations. 

The relation between stellar mass and metallicity in galaxies 
(e.g., Tremonti et al. 2004) suggests that star formation proceeds 
slowly at low metallicity. However, the mass-metallicity rela- 
tion can also be understood as a lower star-formation efficiency 
at low mass, with a relatively more efficient ejection of enriched 
material, or by a metallicity increase through star-formation trig- 
gered by infalling gas. Mannucci et al. (2010) showed that a 
more fundamental relation exists between stellar mass, metal- 
licity, and SFR, where downsizing, outflows, and infalls are the 
major actors. 

Accretion of inflowing cold gas triggers a transient phase in 
which (1) metallicity is diluted by the injection of significant 
amounts of low-metallicity gas; (2) the SFR increases through 
the Schmidt-Kennicutt law; and (3) metallicity increases again 
as star formation proceeds. Although H i companions are often 
observed around star-forming dwarf galaxies (e.g., Taylor 1997; 
Ramya et al. 2009; Pustilnik et al. 2001), their accretion and the 
consequence on star formation is not definitive (e.g., Telles & 
Maddox 2000). Gas inflow is expected to play a more important 
role at high-z. The cosmological inflow from the cosmic web 
could have been so high at high-z that star-formation could not 
catch up in the least massive haloes like IZw 18 (Krumholz & 
Dekel 2012). This metallicity-dependent quenching would result 
in a growth of the Hi envelope, the latter still being observed 
today. Only at z < 2 could a low-mass galaxy like IZw 18 start 
to build up its metallicity. Such a scenario remains compatible 
with the presence of old stars in IZw 18. 

Alternatively (or possibly in parallel), galactic winds could 
expel most of the metals into the IGM through multiple super- 
novae explosions, effectively keeping the ISM metallicity to low 
levels (e.g., Recchi et al. 2004). The most active star-forming 
regions would be responsible for those winds. A major uncer- 
tainty remains however on the fate of the metals observed in 
galactic outflows, since the superbubble expansion might de- 
celerate, and the outflowing material might cool down and re- 
turn to the ISM instead of escaping the gravitational well (e.g., 
Tenorio-Tagle et al. 1999; Mac Low & Ferrara 1999; Rieschick 
2003). Despite the current starburst episode in IZw 18, there is 
no evidence of outflowing hot gas (Sect. 7.1). Furthermore, we 
note that outflows would play a more important role in galaxies 
whose star-formation history is dominated by a series of violent 
bursts, which has probably not been the case in IZw 18 (Legrand 
et al. 2000). 

In summary, although gas flows play a major role in the cos- 
mic evolution of galaxies, their effects in I Zw 1 8 are not yet fully 
constrained and understood. A combination of metal-poor gas 
infall in the early stages of the galaxy, and a chemical downsiz- 
ing could be responsible for the extremely low-metallicity ob- 
served today in the H i and H n regions. Finally, our results do not 
allow putting constraints on the reasons of the current starburst 
episode in IZw 18, which has been ~10-100 times relatively 
stronger over the past ~50 Myr (Aloisi et al. 1999; Legrand et al. 
2000). The sudden onset of star formation could have been trig- 
gered by gas compression from merging events. The disturbed 
H l morphology suggests that several small structures might have 
been interacting (Lelli et al. 2012). Such structures could have 
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been gas-rich dwarf galaxies or sub-DLAs. From the radial ve- 
locity of the neutral metals observed with COS, we did not 
find any evidence of infalling material on the NW region itself; 
the H i gas is, if anything, outflowing from the H n region with 
a velocity of ~20kms _l (see also A03; L04). However, Lelli 
et al. (2012) found a radial inflow/outflow motion on the order 
of ~15kms _1 on larger scales and concluded that interactions 
or mergers are the most likely explanation for the current star- 
burst episode. In general, the FUSE results for other BCDs show 
blue shifts of at most 10-20 km s _l (e.g., Lebouteiller et al. 2009; 
Cannon et al. 2005), confirming that slow radial flows of H i gas 
might be frequent. 

7.3.2. Local conditions for star formation 

We now discuss the star-formation process itself, and the influ- 
ence of metallicity in the ability of clouds to form stars in an 
environment as metal-poor as IZw 18. Our present study puts 
strong constraints on the metallicity of the H t region, which con- 
stitutes the bulk of the mass of the galaxy. The H i region is thus 
key to understand the past star-formation efficiency. 

First, the fact that stars of all masses have formed for several 
Gyrs has some implications for the early metallicity evolution of 
the galaxy. Bromm & Loeb (2003) showed that when the metal- 
licity reaches more than 10 -4 — 10 -3 Z Q , the cooling by the far-IR 
fine-structure lines [Cn] 157yum and [Oi] 63 fim enables cloud 
fragmentation and therefore low-mass star-formation. No dust or 
molecules were considered for this calculation and we note that 
H 2 could provide the most efficient cooling at even lower metal- 
licity (e.g., Omukai et al. 2005). However, the lack of diffuse H 2 
in IZw 18 (Vidal-Madjar et al. 2000) suggests that H 2 cooling is 
not significant in our practical study, except if present in small 
dense clumps where H 2 might be more abundant (see Thuan 
et al. 2005). We thus get the picture of IZw 18 forming stars of 
all masses only a few Gyrs after the Big Bang, when its metal- 
licity reached >10 -4 Z e . The oldest stars in IZw 18 could be the 
remnants of the star-formation episode following this early en- 
richment. The existence of a metallicity threshold for low-mass 
star-formation is however still debated and we note that the pres- 
ence of dust to shield UV photons, and the D/G ratio, ultimately 
related to metallicity via the evolution of the dust-to-metal ratio, 
may be the most important parameters governing the star forma- 
tion process (Schneider et al. 2012; Glover & Clark 2012). 

The importance of the ISM metallicity is fundamental for 
the ability of clouds to form stars, and it may well regulate the 
SFR as long as external gas flows are unimportant. Although the 
observed SFR in BCDs is ~ 10 -2 — 10 M Q yr -1 , this corresponds 
to the current starburst episode, and a much milder SFR could 
have dominated beforehand. For example, in IZw 18 the past 
SFR could have been as low as ~10 -4 M Q yr -1 ; e.g., Legrand 
et al. 2000). It seems at first contradictory with Mannucci et al. 
(2010) who finds that for a given stellar mass, the SFR is anti- 
correlated with metallicity. The fundamental metallicity rela- 
tion found by Mannucci et al. (2010) is, however, biased toward 
emission-line galaxies, and a “pre-starburst” IZw 18 would not 
have been selected by their survey. Furthermore, as discussed in 
the following, it is possible that the extremely low-metallicity 
of this galaxy has a major impact on the star-formation process 
itself. 

The temperature and density of the ISM are essential param- 
eters for star formation. The latter can be inhibited by the in- 
ability of the gas to cool down enough for removing thermal 
pressure and for creating gravitational instabilities. The warm 
temperature in a metal-poor ISM is due to the paucity of efficient 


cooling mechanisms, such as molecules or the far-IR fine- 
structure line [Cn], which dominates the cooling in the diffuse 
atomic clouds. Krumholz (2012) showed that the thermal equi- 
librium timescale is more important than the timescale for con- 
version of atomic to molecular gas in a low-metallicity gas. Star 
formation could proceed before the medium becomes molec- 
ular when the metallicity reaches below a few percents so- 
lar. The metallicity threshold above which thermal equilibrium 
is reached before a few free-fall timescales (once the gas be- 
comes gravitationally unstable) thus depends mostly on the abil- 
ity of the gas to cool down. The Hi region of IZw 18 is an 
ideal testbed for star-formation timescales in a low-metallicity 
atomic gas, since molecular gas has yet to be detected. The 
model in Krumholz (2012) highly depends on the heating mech- 
anisms and their metallicity dependence, and the uncertain heat- 
ing mechanisms in the H 1 region of galaxies such as I Zw 1 8 
(Sect. 5.4) add to the uncertainty on the threshold metallicity 
and thermal equilibrium timescale. Following the method in 
Krumholz (2012), we use our determination of the heating rate 
in the Hi region, and an initial gas density and temperature 
of ~10 cm -3 and ~1000K, respectively. We find that thermal 
equilibrium is reached before a few free-fall timescales only for 
a metallicity greater than ~ 1/100 Z Q . A lower threshold metal- 
licity would be obtained if we consider a significant clumping 
factor. Hence we see that widespread star formation is possi- 
ble in the Hi region of IZw 18, currently holding a metallicity 
of ~ 1/46 Z Q . If such physical conditions prevail in objects as 
metal-poor as < 1/ 1 00 Z Q , the thermal equilibrium timescales are 
rather long and the ability to form stars ultimately depends on the 
survival of the clouds. I Zw 1 8 might have kept a low metallic- 
ity because its ISM remained warm on long timescales. A ma- 
jor uncertainty remains on the presence of small cold clumps 
(cold neutral medium phase) in which star formation would be 
more efficient. The mixing factor of the warm neutral medium 
and cold dense clumps in the FUV absorption lines is difficult 
to quantify. Moreover the ISM pressure in IZw 18 might be too 
low to maintain a multi-phase ISM, so that the neutral gas we 
characterized in this study might be the only thermally-stable 
phase. The temperature from this phase could thus be a crucial 
parameter to understand the star-formation process on large spa- 
tial scales and over cosmic times. 


8. Conclusions 

We present a comprehensive investigation of the H 1 envelope of 
IZw 18. The Hi region holds several fundamental clues to the 
chemical evolution of galaxies. The metallicity of the Hi region 
is set by the successive star-formation episodes over the galaxy 
history, and is modulated by outflows of supernovae products 
and hypothetical inflows of quasi-pristine gas from the cosmic 
web. The abundance of metals in the H 1 region is also a crit- 
ical parameter to understand the present and the future star- 
formation history, as metals provide most of the cooling in the 
H 1 gas reservoir to facilitate and sustain the star-formation pro- 
cess. Our objectives are to characterize the Hi region by (1) its 
physical conditions and (2) its chemical abundances. We com- 
pare the abundances to those derived in the H 11 region (from op- 
tical emission-lines), and discuss the enrichment history of the 
galaxy. 

The Hi envelope of IZw 18 was observed toward the NW 
young stellar cluster in the FUV with the COS instrument 
on HST. The massive stars provide the FUV continuum over 
which absorption-lines from neutral species are superimposed. 
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We modeled the stellar absorption in order to remove its con- 
tribution and isolate the interstellar absorption. We derived the 
column density of H 1 and of heavy elements (C n, C n*, N i, 0 1 , 
Alii, Sin, Pit, Sn, Mnn, Fen, Nin). Models were constructed 
to compute the ionization structure of the nebula. Our models 
predict a large ionized gas volume in the H 11 region, suggesting 
that the H i cavity observed at 21 cm toward NW is due to ioniza- 
tion by the stellar cluster as opposed to blow-out by stellar winds 
and supernovae. We study potential emission components in the 
Hn region that could work out to reduce the apparent optical 
depth of absorption lines. 

The physical conditions in the H i region are drawn from the 
study of the C n* absorption-line. 

- The population of the fine-structure level of C + toward NW 
is compatible with collisions with H° atoms, implying that 
C ii* arises in the H i envelope. 

- We determine observationally that the electron fraction in the 
Hi region is less than ~2%; our models predict ~0.1%. 

- From C n* we derive a cooling rate that is larger than the 
models. An extra heating source is required in the models to 
reconcile them with the observations, leading to an elevated 
electron temperature in the H l region of ~ 1000 K. 

- The cooling rate inferred from Cii* is used to calculate 
the star-formation rate. Our value agrees with other deter- 
minations in the literature, although on the high-end. The 
D/G ratio used for this calculation is a factor of ~2000 times 
lower than the Milky Way value, i.e., a factor of ~60 times 
lower than what is expected if we simply scale D/G with 
metallicity. 

We then infer the dominant ionization stages in the Hi region 
from the models and calculate the chemical abundances in the 
H i region. 

- We find that sulfur is less abundant in the H i region as com- 
pared to the Hn region, by a (log) factor of 0.18 ± 0.1 1. 

- Nitrogen, and possibly carbon, are also underabundant in the 
H i region. 

- The P/S ratio is solar while Al/S is consistent with yields 
from massive stars. The Fe abundance is the same in 
the 2 phases within errors. The [o'/Fe] ratio is consistent with 
production in massive stars. 

- Depletion on dust grains is unlikely to explain the observed 
abundance discontinuity between the Hi and Hn regions. 
This is due to the low dust mass available. 

Finally, we discuss possible scenarios concerning the abundance 
difference in the gas phases, and concerning the enrichment of 
the H i region. 

- Local enrichment in the Hii regions could explain the en- 
hanced N abundance but fails to explain other results, espe- 
cially for S. 

- The presence of relatively lower metallicity outskirts in the 
Hi region is discussed. As much as 50% of the Hi mass 
could be pristine. 

- The metallicity of the Hi region is not null, it is compati- 
ble with the IGM metallicity at z = 1.6-2. 9. IZw 18 could 
be considered as an old DLA with little recent star forma- 
tion. The presence of old stars confirm that IZw 18 is not 
a genuinely young galaxy. The inflow of quasi-pristine gas 
from the cosmic web for a few Gyrs after the Big Bang 
could explain the massive Hi envelope and the quenched 
star-formation in the early stages of the galaxy evolution. 


- The current onset of star formation could be a due the merg- 
ing of dwarfs or sub-DLAs, as suggested by the disrupted 
H i morphology. 

The sample of BCDs whose neutral gas has been investigated 
by FUSE and COS shows that a systematic discontinuity exists 
between the H i and H n region abundances. I Zw 1 8 is the object 
with the lowest discontinuity. Another BCD, SBS 0335-052, de- 
spite having a metallicity similar to I Zw 1 8 does show a signif- 
icant discontinuity which could translate in the H l region being 
as metal-poor as 1/150 Z Q . 

Acknowledgements. V.L. is supported by a CEA/Marie Curie Eurotalents fel- 
lowship. We would like to thank Aurelie Remy for the Herschel dust mass mea- 
surements. Based on observations obtained with the NASA/ESA Hubble Space 
Telescope, which is operated by the Association of Universities for Research in 
Astronomy, Inc., under NASA contract NAS5-26555. S.R.H. and I.H. gratefully 
acknowledge support by the NASA Cosmic Origins Spectrograph program at the 
University of Colorado, Boulder. 


References 

Aloisi, A., Tosi, M„ & Greggio, L. 1999, AJ, 118, 302 
Aloisi, A., Savaglio, S., Heckman, T. M., et al. 2003, ApJ, 595, 760 
Aloisi, A., Annibali, F., Mack, J., et al. 2007, Stellar Populations as Building 
Blocks of Galaxies, 241, 310 

Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481 
Bakes, E. L. O., & Tielens, A. G. G. M. 1994, ApJ, 427, 822 
Barbuy, B. 1988, A&A, 191. 121 

Bhattal, A. S„ Francis, N„ Watkins, S. J., & Whitworth, A. P. 1998, MNRAS, 
297, 435 

Bowen, D. V., Jenkins, E. B., Pettini, M., & Tripp, T. M. 2005, ApJ, 635, 880 
Brinchmann, J., Kunth, D., & Durret, F. 2008, A&A, 485, 657 
Bromm, V., & Loeb, A. 2003, Nature, 425, 812 
Brooks, A. M., Governato, F., Booth, C. M., et al. 2007, ApJ, 655, L17 
Brown, T. M., Heap, S. R., Hubeny, I., Lanz, T., & Lindler, D. 2002, ApJ, 579, 
L75 

Cannon, J. M., Skillman, E. D.. Sembach, K. R., & Bomans, D. J. 2005, ApJ, 
618, 247 

Castelli, F„ & Kurucz, R. L. 2003, ApJ, 210 
Cen, R„ & Ostriker, J. P 1999, ApJ, 519, L109 

Contreras Ramos, R., Annibali, F., Fiorentino, G., et al. 201 1, ApJ, 739, 74 
Crawford, M. K„ Genzel, R„ Townes, C. H„ & Watson, D. M. 1985, ApJ, 291. 
755 

Dufour, R. J., & Hester, J. J. 1990, ApJ. 350, 149 
Elmegreen, B. G. 201 1, ApJ, 51, 59 

Federman, S. R., Brown, M.. Torok, S., et al. 2007, ApJ, 660, 919 
Ferland, G. J., Korista, K. T„ Vemer, D. A., et al. 1998, PASP. 1 10, 761 
Fukuda, N„ & Hanawa, T. 2000, ApJ, 533, 911 

Galliano, F., Madden, S. C., Jones, A. P, Wilson, C. D., & Bernard, J. P. 2005, 
A&A, 434, 867 

Galliano, F„ Dwek, E„ & Chanial, P. 2008, ApJ, 672, 214 
Garnett, D. R., Skillman, E. D., Dufour, R. J., et al. 1995, ApJ, 443, 64 
Garnett, D. R„ Skillman, E. D„ Dufour, R. J., & Shields, G. A. 1997, ApJ, 481. 
174 

Giovanelli, R„ Haynes, M. P, Kent, B. R„ et al. 2005, AJ, 130, 2598 
Glover, S. C. O., & Clark, P. C. 2012, MNRAS, 421, 1 16 
Goldsmith, P. F., Langer, W. D., Pineda, J. L., & Velusamy, T. 2012, ApJS, 203, 
13 

Green, J. C., Froning, C. S., Osterman, S., et al. 2011, ApJ, 744, 60 
Grimes, J. P., Heckman, T.. Aloisi, A., et al. 2009, ApJ, 181, 272 
Gry, C., Lequeux, J., & Boulanger, F. 1992, A&A, 266, 457 
Habing, H. J. 1968, ApJ, 19, 421 

Heap, S., Malumuth, E., Landsman, W., & Lindler, D. 2013, ApJ, submitted 
Hebrard, G., Lemoine, M., Vidal-Madjar, A., et al. 2002, ApJ, 140, 103 
Heckman, T. M., Sembach, K. R., Meurer, G. R., et al. 2001, ApJ, 554, 1021 
Heckman, T. M., Norman, C. A., Strickland, D. K., & Sembach, K. R. 2002, 
ApJ, 577, 691 

Henry, R. B. C., Edmunds, M. G„ & Koppen, J. 2000, ApJ, 541, 660 
Henry, R. B. C„ Nava, A., & Prochaska, J. X. 2006, ApJ, 647, 984 
Howk, J. C., Wolfe, A. M„ & Prochaska, J. X. 2005, ApJ, 622, L8 1 
Hubeny, I., & Lanz, T. 1995, ApJ, 439, 875 

Hubeny, I., Heap, S. R„ & Lanz, T. M. 2000, AAS 197th meeting, 32, 718 
Hulsbosch, A. N. M„ & Wakker, B. P. 1988, ApJ, 75, 191 
Hunter, D. A., & Thronson, H. A. J. 1995, ApJ, 452, 238 


A16, page 24 of 32 


V. Lebouteiller et al.: Chemical enrichment and physical conditions in IZw 18 


Hunter, I., Brott, I., Lennon, D. J., et al. 2008, ApJ, 676, L29 

Izotov, Y. I., & Thuan, T. X. 1999, ApJ, 51 1, 639 

Izotov, Y. I., & Thuan, T. X. 2007, ApJ, 665, 1115 

Izotov, Y. I., Thuan, T. X., & Lipovetsky, V. A. 1997, ApJ, 108, 1 

Izotov, Y. I., Chaffee, F. H„ Foltz, C. B„ et al. 1999, ApJ, 527, 757 

Izotov, Y. I., Chaffee, F. H„ Foltz, C. B„ et al. 2001a, ApJ, 560, 222 

Izotov, Y. I., Schaerer, D„ & Charbonnel, C. 2001b, ApJ, 549, 878 

Izotov, Y. I., Thuan, T. X., & Guseva, N. G. 2005, ApJ, 632, 210 

Jenkins, E. B. 1971, ApJ. 169, 25 

Jenkins, E. B. 2009, ApJ, 700, 1299 

Jenkins, E. B. 2013, ApJ. 764, 29 

Jenkins, E. B„ Gry, C., & Dupin, O. 2000a, A&A, 354, 253 

Jenkins, E. B„ Oegerle, W. R„ Gry, C., et al. 2000b, ApJ, 538, L81 

Jonsell, K., Edvardsson, B., Gustafsson, B., et al. 2005, ApJ, 440, 321 

Kennicutt, Jr, R. C. 1998, ApJ, 498, 541 

Kimura, H., Mann, I., & Jessberger, E. K. 2003, ApJ, 582, 846 

Kobulnicky, H. A., & Skillman, E. D. 1997, ApJ, 489, 636 

Krumholz, M. R. 2012, ApJ, 759, 9 

Krumholz, M. R„ & Dekel, A. 2012, ApJ, 753, 16 

Kunth, D„ & Ostlin, G. 2000, A&ARv, 10, 1 

Kunth, D„ & Sargent, W. L. W. 1986, ApJ, 300, 496 

Kunth, D„ Lequeux, J„ Sargent, W. L. W„ & Viallefond, F. 1994, A&A, 282, 
709 

Lanz, T„ & Hubeny, I. 2003, ApJ, 146, 417 
Lanz, T„ & Hubeny, I. 2007, ApJ. 169, 83 
Lebouteiller, V. 2005, Ph.D. Thesis 

Lebouteiller, V., Kunth, D.. Lequeux, J., et al. 2004, A&A, 415, 55 
Lebouteiller, V., Kuassivi, & Ferlet, R. 2005, A&A, 443, 509 
Lebouteiller, V., Kunth, D., Lequeux, J.. et al. 2006, A&A, 459, 161 
Lebouteiller, V., Bernard-Salas, J., Brandi, B., et al. 2008, ApJ, 680, 398 
Lebouteiller, V., Kunth, D„ Thuan, T. X., & Desert, J.-M. 2009, A&A, 494, 915 
Lecavelier des Etangs, A., Desert, J.-M., Kunth, D., et al. 2004, A&A, 413, 131 
Lee, H. T„ & Chen, W. P. 2007, ApJ, 657, 884 
Lee, H„ & Skillman, E. D. 2004, ApJ, 614, 698 

Legrand, F., Kunth, D., Roy, J.-R., Mas-Hesse, J. M., & Walsh, J. R. 2000, A&A, 
355, 891 

Lehner, N„ Wakker, B. P„ & Savage, B. D. 2004, ApJ, 615, 767 
Lejeune, T., & Schaerer, D. 2001, ApJ, 366, 538 

Lelli, F., Verheijen, M., Fraternali, F., & Sancisi, R. 2012, A&A, 537, A72 
Lemoine, M., Vidal-Madjar, A., Hebrard, G., et al. 2002, ApJ, 140, 67 
Lequeux, J., & Viallefond, F. 1980, A&A, 91, 269 
Lodders, K. 2003, ApJ, 591, 1220 

Lopez Sanchez, A. R., Esteban, C., Garcia Rojas. J., Peimbert, M., & Rodriguez, 
M. 2007, ApJ, 656, 168 

Mac Low, M. M., & Ferrara, A. 1999, ApJ, 513, 142 
Maeder, A., & Meynet, G. 2012, ApJ, 84, 25 
Madden, S. C., Geis, N„ Genzel, R„ et al. 1993, ApJ, 407, 579 
Mannucci, F., Cresci, G., Maiolino, R., Marconi, A., & Gnerucci, A. 2010, 
MNRAS, 408, 2115 

Meynet, G., & Maeder, A. 2002a, A&A, 390, 561 

Meynet, G„ & Maeder, A. 2002b, A&A, 381, L25 

Moos, H. W„ Cash, W. C„ Cowie, L. L„ et al. 2000, ApJ, 538, LI 

Morton, D. C. 2003, ApJ, 149, 205 

Nava, A., Casebeer, D., Henry, R. B. C., & Jevremovic, D. 2006, ApJ, 645, 1076 
Noeske, K. G„ Guseva, N. G„ Fricke, K. J., et al. 2000, A&A, 361, 33 
Omukai, K., Tsuribe, T., Schneider, R., & Ferrara, A. 2005, ApJ, 626, 627 
Papaderos, P, & Ostlin, G. 2012, A&A, 537, A126 
Peimbert, A., & Peimbert, M. 2010, ApJ, 724, 791 
Pequignot, D„ 2008, A&A, 478, 371 

Petrosian, A. R., Boulesteix, J., Comte, G., Kunth, D., & Lecoarer, E. 1997, 
A&A, 318, 390 


Pettini, M„ & Lipman, K. 1995, A&A, 297, L63 

Pettini. M., Steidel, C. C., Adelberger, K. L., Dickinson, M., & Giavalisco, M. 
2000, ApJ, 528, 96 

Podobedova, L. I., Kelleher, D. E., & Wiese, W. L. 2009. J. Phys. Chem. Ref. 
Data, 38, 171 

Pottasch, S. R., Wesselius, P. R., & van Duinen, R. J. 1979, A&A, 74, L15 
Prochaska, J. X., & Wolfe, A. M. 1999, ApJ, 121, 369 

Pustilnik, S. A., Brinks, E., Thuan, T. X., Lipovetsky, V. A., & Izotov, Y. I. 2001, 
AJ, 121, 1413 

Rafelski, M., Wolfe, A. M., Prochaska, J. X., Neeleman, M., & Mendez, A. J. 
2012, ApJ, 755, 89 

Ramya, S., Kantharia, N. G„ & Prabhu, T. P. 2009, ASP Conf. Set., 407, 1 14 
Rana, N. C. 1991, ARA&A, 29, 129 

Recchi, S., Matteucci, F., D’Ercole, A., & Tosi, M. 2004, A&A, 426, 37 
Rieschick, A. 2003, ApJ, 284, 861 
Sargent, W. L. W„ & Searle, L. 1970, ApJ, 162, L155 
Savage, B. D„ & Sembach, K. R. 1996, ARA&A, 34, 279 
Schlegel, D. J., Finkbeiner, D. P, & Davis, M. 1998, ApJ, 500, 525 
Schneider, N., Glisten, R., Tremblin, P., et al. 2012, A&A, 542, L18 
Schulte Ladbeck, R. E., Rao, S. M., Drozdovsky, 1. O., et al. 2004, ApJ, 600, 613 
Sembach, K. R„ Howk, J. C., Ryans, R. S. I., & Keenan, F. P. 2000, ApJ, 528, 
310 

Shapley, A. E„ Steidel, C. C„ Pettini, M„ & Adelberger, K. L. 2003, ApJ, 588, 
65 

Skillman, E. D„ & Kennicutt, R. C. J. 1993, ApJ, 411, 655 
Sofia, U. J., & Jenkins, E. B. 1998, ApJ, 499, 951 

Sofia, U. J., Cardelli, J. A., Guerin, K. P„ & Meyer, D. M. 1997, ApJ, 482, L105 
Sofia, U. J., Fitzpatrick, E. L„ & Meyer, D. M. 1998, ApJ, 504, L47 
Taylor, C. L. 1997, ApJ, 480, 524 

Telfer, R. C.. Kriss, G. A., Zheng, W., Davidsen, A. F., & Tytler, D. 2002, ApJ, 
579, 500 

Telles, E„ & Maddox, S. 2000, MNRAS, 311, 307 

Tenorio-Tagle, G., Silich, S. A., Kunth, D., Terlevich, E., & Terlevich, R. 1999, 
MNRAS, 309, 332 

Thuan, T. X., & Izotov, Y. I. 2005, ApJ, 161, 240 
Thuan, T. X., & Martin, G. E. 1981, ApJ, 247, 823 
Thuan, T. X., Izotov, Y. I., & Lipovetsky, V. A. 1996, ApJ. 463, 120 
Thuan, T. X., Izotov, Y. I., & Lipovetsky, V. A. 1997, ApJ, 477, 661 
Thuan, T. X., Lecavelier Des Etangs, A. L. d., & Izotov, Y. I. 2002, ApJ, 565, 
941 

Thuan, T. X., Bauer, F. E., Papaderos, P., & Izotov, Y. I. 2004, ApJ, 606, 213 
Thuan, T. X., Lecavelier des Etangs, A., & Izotov, Y. I. 2005, ApJ, 621, 269 
Tremonti, C. A., Heckman, T. M., Kauffmann, G., et al. 2004, ApJ, 613, 898 
van Zee, L„ Westpfahl, D„ & Haynes, M. P. 1998, AJ, 1 15, 1000 
Vanhala, H. A. T„ & Cameron, A. G. W. 1998, ApJ, 508, 291 
Verner, D. A., Ferland, G. J., Korista, K. T„ & Yakovlev, D. G. 1996, ApJ, 465, 
487 

Vidal-Madjar, A., Kunth, D., Lecavelier Des Etangs, A. L. d., et al. 2000, ApJ, 
538, L77 

Wakker, B. P„ & van Woerden, H. 1997, ARA&A, 35, 217 
Walsh, J. R„ & Roy, J.-R. 1989, MNRAS, 239, 297 
Walsh, J. R„ & Roy, J.-R. 1993, MNRAS, 262, 27 
Weingartner, J. C., & Draine, B. T. 2001, ApJ, 134, 263 
Wolfe, A. M„ Gawiser, E„ & Prochaska, J. X. 2003a, ApJ, 593, 235 
Wolfe, A. M„ Prochaska, J. X., & Gawiser, E. 2003b, ApJ, 593, 215 
Wolfire, M. G., Hollenbach, D„ McKee, C. F„ Tielens, A. G. G. M„ & Bakes, E. 
L. O. 1995, ApJ, 443, 152 

Wood, B. E„ Linsky, J. L„ Hebrard, G., et al. 2002, ApJ, 140, 91 
Woosley, S. E„ & Weaver, T. A. 1995, ApJ, 101, 181 
Wu, Y„ Charmandaris, V., Hao, L.. et al. 2006, ApJ, 639, 157 
Zwicky, F. 1966, ApJ, 143, 192 


Pages 26 to 32 are available in the electronic edition of the journal at http : //www . aanda . org 


A16, page 25 of 32 


A&A 553, A16 (2013) 


Appendix A: Line fitting 


1.2 

1.0 

0.8 

0.6 

0.4 

0.2 

0.0 


CO_ 

o o a5 



Wavelength [A] 


Wavelength [A] 


CVJ, — cm cm 

i_ _ *_ 

CL CL CL O 



-■cm co cmcm ' — 1 - t— t- cm 


CM CM CO CM -I- 1- i- - 


LLJ,,V UJ 


ooo 

u; 




Wavelength [A] 


Fig.A.l. Line profile fitting for the 1 1 523 (top), 12 028 {middle), and 1 1 579 datasets (bottom). For display purposes, only the strongest lines 
from IZw 18 (component 1, also indicated by the vertical dotted lines), from the Milky Way (component 2), and from the high-velocity cloud 
(component 3) are labeled. The thin green line shows the individual line profiles. The solid thick red line shows the fit to the data while the dotted 
red line shows the line profiles assuming the FI i region abundances are equal to the H n region abundances. 
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Fig.A.2. See Fig. A.l for the plot description. The unlabelled line on the far right is Sim /11206.5 arising in IZw 18. 
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Fig.A.3. See Fig. A.l for the plot description. 


A16, page 28 of 32 


V. Lebouteiller et al.: Chemical enrichment and physical conditions in IZw 18 



ooooo oooooozoo 


U4,,u uyyy igg 




OOOOO oooooozoo 



1310 1312 1314 1316 1318 1320 1322 1324 1325 1330 1335 1340 1345 1350 

Wavelength [A] Wavelength [A] 



CO oooo o ooooo zoo 


i , y tgy 



Fig.A.4. See Fig. A.l for the plot description. 
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Fig.A.5. See Fig. A.l for the plot description. 
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Fig.A.6. See Fig. A.l for the plot description. 
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Fig.A.7. See Fig. A.l for the plot description. 
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